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Our  increasing  dependence  on  electricity  brings  awareness  of  the 
value  of  reliable  service  and  the  inconvenience  and  loss  to  consumers 
that  may  stem  from  interruptions  to  supply  and  load  shedding.  Since  the 
cost  of  reliability  is  significant,  decision-makers  need  to  choose 
prices  and  capacity  so  as  to  balance  additional  benefits  against  the 
additional  cost  of  reliability.  Instead  of  trying  to  capture  both  the 
benefit  and  cost  sides  of  alternative  reliability  levels,  this  research 
will  concentrate  on  the  understanding  of  the  cost  side  and  the  implica- 
tions of  alternative  reliability  indices. 

This  study  differs  from  previous  economic  studies  in  three 
respects.  The  reliability  of  power  supply  becomes  an  essential  variable 


in  the  analysis  of  electricity  costs.  Two  reliability  indices  widely 
used  by  engineers,  loss  of  load  probability,  LOLP,  and  expected  demand 
not  served,  e(DNS),  are  considered  here.  The  second  contribution  is 
that  demand  and  supply  are  both  treated  as  stochastic  variables.  The 
third  aspect  of  the  study  warranting  attention  is  the  use  of  pseudo 
data.  The  numbers  allow  us  to  apply  econometric  estimation  techniques 
to  determine  the  cost  of  reliability. 

Each  of  the  indices  has  its  own  strengths  and  weaknesses,  and 
cannot  individually  provide  a complete  description  of  outage  patterns, 
nor  of  the  value  of  lost  consumption.  The  estimated  indices  are  only 
significant  in  relative,  rather  than  absolute  terms. 

Due  to  lack  of  sufficient  and  reliable  data,  pseudo  data  based  on 
underlying  engineering  and  economic  considerations  were  used  in  this 
study.  The  model  we  propose  is  a set  of  two  equations  which  consist  of 
a reliability  equation  and  a cost  equation.  A new  concept,  reliability 
elasticity  of  cost,  was  introduced  to  measure  the  sensitivity  between 
cost  and  reliability.  The  statistical  testing  of  estimated  coefficients 
suggests  that  the  e(DNS)  model  provides  the  greatest  explanatory  power. 

As  to  future  research,  we  recommend  three  promising  areas.  They 
are  the  benefit  side  of  reliability,  interruptibl e service  pricing,  and 
generating  unit  size  selection. 
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CHAPTER  1 
INTRODUCTION 

Electricity  is  a dominant  feature  of  industrialized  societies, 
serving  as  an  input  to  industrial  and  commercial  production  and  the 
major  source  of  energy  for  the  household  sector.  Consumption  of  elec- 
tricity has  consistently  increased  over  the  last  few  decades,  doubling 
at  least  every  ten  years  in  developed  countries.  Utility  managers  are 
beginning  base  load  construction  programs  which  will  increase  capacity 
by  35  percent  through  the  late  1980s.  Since  generating  plants  take  a 
decade  or  more  to  plan  and  build,  proper  forecasting  and  the  selection 
of  appropriate  reliability  goals  are  key  tasks  facing  utilities  and 
regulators.  Our  increasing  dependence  on  electricity  brings  awareness 
of  the  value  of  reliable  service  and  the  inconvenience  and  loss  to 
consumers  due  to  supply  interruptions  and  load  shedding.  Yet  we  know 
that  the  cost  of  reliability  is  significant,  requiring  decision  makers 
to  choose  prices  and  capacity  to  balance  additional  benefits  against 
additional  costs  of  reliability. 

1 . 1 Purpose  of  Study 

According  to  Electrical  Power  Research  Institute  (EPRI)  studies 
(1978),  almost  one  third  of  investment  in  power  systems  is  utilized  in 
standby  plant  and  redundant  transmission  networks  to  achieve  a higher 
reliability  of  supply.  Due  to  the  capital  intensive  nature  of  power 
systems,  evaluation  of  reliability  economics  is  of  primary  importance  in 
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importance  in  achieving  optimum  allocation  of  resources  within  the 
electricity  supply  industry.  Furthermore,  the  issue  of  reliability  had 
arisen  out  of  a concern  on  the  part  of  the  public  and  of  regulators  that 
utilities  may  be  incurring  excessive  costs  for  reliability.  This  con- 
cern has  come  at  a time  when  utilities  are  seeking  to  curb  expenditures, 
and  we  might  well  be  concerned  that  utilities  are  cutting  so  deeply  into 
their  construction  programs  that  over  the  long-run  the  reliability  of 
electric  service  will  be  severely  jeopardized. 

From  the  engineering  standpoint,  system  reliability  is  one  of  the 
major  factors  in  the  planning,  design,  and  operation  of  power  systems, 
since  reliability  criteria  have  a major  impact  on  the  requirement  of  new 
facilities  and  system  operating  security.  The  criteria  ultimately 
affect  the  quality  of  service  and  the  rates  customers  pay.  Utilities 
are  therefore  under  pressure  to  make  more  precise  evaluations  of  relia- 
bility criteria  and  the  associated  costs  and  benefits. 

1.2  Outline  of  Reliability  Issues 

Power  system  reliability  is  one  aspect  of  the  wider  topic  of  the 
effectiveness  of  energy  supply  systems  as  a whole.  Some  of  the  inter- 
related issues  include 

1.  Concepts  of  reliability 

a.  Behavioral  impl ications--risk  indices  in  planning  and 
operation 

b.  Technical  analysis--methodologies  for  evaluating 
reliability 

c.  Scheduling  algorithms  for  maintenance  and  repair 

2.  Load  forecasting  and  demand  uncertainties 
Uncertainties  in  primary  energy  availability 
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4.  Load  management  and  reliability 

a.  Peak  load  or  time-of-use  pricing 

b.  Interruptible  service  pricing 

c.  Direct  load  control 

5.  Determination  of  optimum  reliability 

a.  Data  collection  and  preparation  for  decision-makers 

b.  Cost  of  reliability  and  its  estimation 

c.  Worth  (benefit)  of  reliability  and  its  estimation 

6.  Environmental  and  socio-economic  consequences  of  non-price 
rationing  (such  as  blackouts) 

a.  Defining  and  quantifying  the  kinds  of  environmental, 
economic,  and  social  effects 

b.  Determining  what  measures  might  be  taken  by  utility 
companies  and  their  customers  during  electricity  short- 
ages to  avoid  costs 

7.  Costs  and  benefits  of  over  and  under  capacity  in  the  genera- 
tion of  electric  power  (Are  the  costs  of  being  underprepared 
for  future  electric  demands  significantly  greater  than  the 
costs  of  being  overprepared?) 

1.3  The  Focus  of  This  Study 

Cost-benefit  analysis  and  its  employment  in  the  planning  of  the 
power  system  has  not  advanced  in  the  way  it  has  in  other  public  services 
(e.g.,  transportation),  partly  due  to  the  great  difficulty  in  obtaining 
comprehensive  assessments  and  measures  of  reliability  at  the  generation 
system  level,  down  to  outage  due  to  distribution  system  failure.  Still 
more  difficult  is  the  assessment  and  inconvenience  of  lost  consumer 
valuations  due  to  interruptions.  It  is  possible,  however,  to  assess  the 
contribution  of  different  schemes  toward  the  continuity  of  supply  and  to 
choose  the  one  with  the  maximum  contribution  within  "reasonable  cost." 
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The  problem  still  remains  as  to  determining  the  reasonable  cost  and  of 
choosing  a reliability  index  that  is  conceptually  sound. 

Systems  planners  do  not  aim  at  a completely  uninterruptible  supply, 
although  the  constraint  or  reliability  criterion  is  usually  taken  as  a 
given  rather  than  analyzed  as  a determinant  of  revenue  requirements. 
Ideally,  we  search  for  a point  of  "optimum  reliability";  that  is,  point 
E in  Figure  1-1.  Two  factors  must  be  balanced:  the  reliability  mar- 
ginal cost  curve  to  the  electric  utility  industry  and  the  marginal 
benefit  of  reliability.  The  latter  concept  depends  on  the  marginal  cost 
of  interruptions  to  the  consumer  (marginal  value  of  the  demand  not 
served).  As  construction  and  financial  costs  have  increased,  measure- 
ment of  the  costs  and  benefits  of  reliability  have  become  increasingly 
important.  Estimation  of  the  cost  curve  is  conceptually  straightforward 
and  shall  be  the  focus  here.  Even  this  one-sided  approach  requires  the 
resolution  of  some  extremely  complex  issues. 

As  to  the  marginal  benefit  of  reliability,  one  analytical  approach 
involves  measuring  the  resulting  reduction  in  consumer  surplus  loss: 
outage  costs  are  reflected  in  the  effects  on  the  production  of  final 
goods  and  services  in  various  sectors.  Conceptually,  a marginal  benefit 
curve  of  reliability  can  be  constructed,  but  few  economists  have  ad- 
dressed this  issue. 

Instead  of  trying  to  capture  both  benefit  and  costs  sides  of  alter- 
native reliability  levels,  this  research  will  sharpen  our  understanding 
of  the  cost  of  reliability  and  implications  of  alternative  reliability 
indices.  In  Chapter  3,  we  first  discuss  the  meanings  and  implications 
of  alternative  reliability  indices,  including  loss  of  load  probability 
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Figure  1-1.  Determination  of  Optimal  Reliability 
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and  expected  demand  not  served.  The  cost  characteristics  of  the  utility 
industry  are  then  reviewed  to  facilitate  the  development  of  a highly 
stylized  characterization  of  these  costs. 

The  framework  is  extended  in  the  next  chapter  to  simulate  the 
electricity  systems.  The  simulations  allow  the  estimation  of  cost  of 
reliability  equations.  In  Chapter  5,  a cost  model  for  the  estimation  of 
reliability  costs  is  presented.  The  econometric  results  are  then  pre- 
sented and  discussed. 

1.4  Summary  of  Conclusions 

Since  the  cost  of  reliability  is  significant,  decision  makers  need 
to  choose  prices  and  capacity  so  as  to  balance  additional  benefits 
against  the  additional  cost  of  reliability.  Instead  of  trying  to  cap- 
ture both  the  benefit  and  cost  sides  of  alternative  reliability  levels, 
this  research  will  concentrate  on  understanding  the  cost  side  and  the 
implications  of  alternative  reliability  indices. 

This  study  differs  from  previous  economic  studies  in  three 
respects: 

1)  The  reliability  of  power  supply  is  an  essential  variable  in 
the  analysis  of  electricity  costs.  Two  reliability  indices 
widely  used  by  engineers,  LOLP  and  e(DNS),  are  considered 
here. 

2)  Demand  and  supply  are  both  treated  as  stochastic  variables. 

3)  Pseudo  data  are  used  to  apply  econometric  estimation  tech- 
niques to  determine  the  cost  of  reliability. 

Each  of  the  indices  has  its  own  strengths  and  weaknesses  and  cannot 
individually  provide  a complete  description  of  outage  pattern  nor  of  the 
value  of  lost  consumption.  The  estimated  indices  are  only  significant 
in  relative,  rather  than  absolute  terms.  Since  the  indices  are 
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generally  neither  comparable  nor  consistent,  they  yield  different  rank- 
ings of  systems  by  different  reliability  indices. 

Due  to  lack  of  sufficient  and  reliable  data,  pseudo  data  based  on 
underlying  engineering  and  economic  considerations  were  used  in  this 
study.  A simulation  model  to  generate  such  data  is  presented  in  Chapter 
4.  We  then  use  alternative  regression  models  to  estimate  the  cost  of 
reliability  in  Chapter  5. 

The  model  we  proposed  is  a set  of  two  equations  which  consist  of  a 
reliability  equation  and  a cost  equation,  appearing  in  a recursive 
system.  A new  concept,  reliability  elasticity  of  cost,  was  introduced 
to  measure  the  sensitivity  of  cost  to  different  reliability  levels.  The 
statistical  testing  of  estimated  coefficients  suggests  that  the  e(DNS) 
model  provides  the  best  estimate. 

As  to  future  research,  we  recommend  three  promising  areas:  the 
benefit  side  of  reliability,  interruptible  service  pricing,  and  gener- 
ating unit  size  selection. 


CHAPTER  2 

REVIEW  OF  THE  LITERATURE 


One  of  the  difficult  problems  faced  by  electric  supply  system 
planners  is  that  of  deciding  how  far  they  are  justified  in  increasing 
the  investment  in  their  equipment  to  improve  service  reliability.  Both 
engineers  and  economists  are  interested  in  this  important  issue. 
Naturally,  engineering  studies  are  engineering-technical  efficiency 
oriented.  This  line  of  research  is  best  represented  by  the  work  of  Roy 
B i 1 1 i ng ton  and  Robert  Sullivan.  Economic  studies,  in  contrast,  are 
economic  efficiency-oriented  and  use  the  concept  of  economic  value. 
Illustrative  studies  are  those  by  Brown  and  Johnson,  Crew  and 
K1 ei ndorfer , and  Munasinghe  and  Gellerson.  In  these  economic  studies, 
there  are  two  different  approaches  to  determining  an  optimal  reliability 
level.  A detailed  discussion  of  these  approaches  will  be  presented  in 
this  chapter. 

The  analysis  of  peak  load  pricing  under  the  deterministic  case, 
which  has  now  become  standard,  produces  the  rule  that  prices  should 
equal  marginal  opportunity  costs  or  the  price  necessary  to  restrict 
demand  to  capacity.  The  first  study  to  introduce  elements  of  risk  into 
the  analysis  was  by  Brown  and  Johnson  (1969).  In  their  model,  the 
social  welfare  function  that  is  maximized  is  the  algebraic  difference 
between  expected  willingness  to  pay  and  expected  cost.  They  concluded 
that  the  optimal  price  will  always  be  lower  and,  with  linear  demand,  the 
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optimal  output  will  generally  be  higher  than  their  counterparts  in  the 
riskless  model  of  traditional  theory,  regardless  of  the  manner  in  which 
the  disturbance  term  enters  the  analysis. 

2.1  Approaches  to  Reliability 

Generally,  there  are  two  different  approaches  to  determining  the 
optimal  reliability  level,  the  welfare-maximum  approach  and  the 
Mu  nasi nghe/Gel 1 erson  benefit/cost  (B/C)  approach.  We  can  summarize  the 
complex  analysis  of  the  welfare-maximum  approach  in  a very  simplified 
form  as  follows:  Economic  efficiency  would  be  maximized  if  price  is  set 
equal  to  marginal  cost,  determining  the  capacity  level  and  reliability 
level  in  the  process.  The  B/C  approach  optimizes  the  reliability  level 
by  comparing  the  social  benefits  and  costs  of  change  in  power  system 
reliability.  These  costs  (supply  side)  can  be  determined  from  engi- 
neering considerations,  such  as  technical  constraints  regarding  capacity 
availability.  On  the  demand  side,  the  benefits  to  electricity  users 
consist  of  cost  savings  from  averted  power  failure,  which  may  be  mea- 
sured by  the  disruption  of  the  output  streams  owing  to  idle  input  fac- 
tors and  spoilage. 

As  can  be  seen  in  Table  2-1,  there  are  two  major  differences  be- 
tween the  two  approaches.  First,  in  the  welfare-maximum  approach,  elec- 
tricity is  treated  as  a good  which  directly  provides  consumers  with 
satisfaction.  Therefore,  outage  costs  are  measured  primarily  in  terms 
of  the  resulting  reduction  in  consumer  surplus  (plus  rationing  costs  in 
the  case  of  Crew  and  K1  eindorfer) . However,  in  the  Munasinghe  and 
Gel  1 erson  B/C  approach,  electricity  is  viewed  as  being  an  intermediate 
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product  which  is  used  to  produce  final  goods  that  are  demanded  by  con- 
sumers. Outage  costs  are  measured  in  terms  of  the  effects  on  the  pro- 
duction of  final  goods  and  services  in  various  sectors  of  the  economy. 
These  approaches  are  not  mutually  exclusive;  each  can  throw  light  on 
this  complex  topic. 

A second  difference  between  the  approaches,  given  a rationing  or 
load  shedding  scheme,  is  that  the  welfare  approach  maximizes  net  welfare 
by  satisfying  joint  first  order  conditions  with  respect  to  prices  and 
capacity.  Partial  approaches  which  attempt  to  optimize  reliability  in 
the  presence  of  given  prices  cannot  ensure  that  resulting  net  benefits 
are  minimized.  However,  in  practice,  electricity  tariffs  are  often  not 
readily  subject  to  change  in  the  short  run.  Therefore,  M&G  developed  a 
method  of  optimizing  reliability  by  weighing  the  costs  and  benefits 
associated  with  alternative  levels  of  reliability,  given  fixed  prices. 
Once  the  optimal  reliability  level  is  determined  in  the  first  round, 
tariffs  can  be  revised  to  reflect  any  changes  in  the  cost  of  supplying 
electricity  implied  by  the  new  reliability  level.  With  this  new  level 
of  tariffs  and  the  resulting  revised  demand  forecast,  reliability  can  be 
re-optimized.  The  process  can  be  continued  until  it  determines  the 
internally  consistent  combination  of  tariffs  and  reliability  that  maxi- 
mizes the  net  social  benefits  of  supplying  electricity. 

2.2  Studies  of  Reliability 

Much  of  the  literature  of  public  enterprise  pricing  ignores  sto- 
chastic phenomena  affecting  demand  and  available  capacity.  For  example, 
the  problem  of  peak  load  pricing  for  public  utilities  has  been  thor- 
oughly explored  under  conditions  of  certainty.  As  mentioned  at  the 
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beginning  of  this  chapter,  the  analysis  of  the  deterministic  case,  which 
has  now  become  standard,  produces  the  rule  that  price  should  equal 
whichever  is  the  higher,  either  marginal  running  costs  or  the  price 
necessary  to  restrict  demand  to  capacity.  A paper  by  Brown  and  Johnson 
(1969),  hereinafter  referred  to  as  BJ,  first  introduced  elements  of  risk 
into  the  analysis. 

BJ  extended  the  problem  to  situations  where  demand  is  subject  to  a 
random  disturbance.  It  was  assumed  that  a producer  is  directed  to  maxi- 
mize expected  social  welfare  and  must  announce  a price  and  capacity 
output  before  actual  demand  is  known.  Leland  (1972)  suggested  that  the 
electric  power  industry  is  a good  example  of  price  setting,  where  prices 
are  set  before  demand  (D)  is  known  and  then  output  is  adjusted  to  meet 
demand.  However,  when  demand  exceeds  capacity,  no  further  output  can  be 
produced.  In  the  short  run,  the  utility  cannot  adjust  capacity  to  meet 
demand  since  capacity  adjustments  require  too  much  time.  This  point 
underscores  the  issues  of  electricity  reliability  which  can  arise. 

In  the  BJ  model,  the  arguments  of  the  welfare  function  to  be  maxi- 
mized include  the  expected  value  of  willingness  to  pay,  the  expected 
value  of  variable  cost,  and  the  capacity  cost  known  with  certainty.  We 
can  define  the  welfare  function  as  follows: 

E(W)  = E(wi 1 1 i ngness  to  pay)  - E(average  variable  cost  • sales) 

- E(capacity  cost) 

= E(consumers 1 surplus  + total  revenue  - total  cost).  (2-1) 

Consumers'  surplus  refers  to  the  usual  consumers'  surplus  triangle 
under  the  demand  curve.  As  can  be  seen  in  Figure  2-1,  when  quantity 
demanded  (given  PQ)  exceeds  capacity  (Z),  part  of  the  triangle  is  not 
attained.  Any  portion  of  the  triangle  beyond  the  capacity  level  is 
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Figure  2-1.  Lost  Consumer  Surplus 
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lost,  since  demand  in  this  region  is  not  satisfied.  This  portion, 
denoted  by  L,  represents  lost  consumer  surplus.  Here  Z represents 
capacity  and  PQ  is  the  going  rate.  This  definition  of  L requires  that 
consumers  are  costlessly  ranked  according  to  their  willingness  to  pay. 
When  demand  exceeds  capacity,  those  consumers  with  the  highest  valuation 
are  served  first.  In  practice,  ranking  consumers  is  not  a costless 

activity,  so  the  lost  consumer  surplus  when  demand  is  not  served  may  be 
greater  than  L. 

Figure  2-2  illustrates  two  other  ways  to  ration  the  excess  demand, 
where  q*  is  energy  demanded.  The  demand  not  served  is  q*  - Z.  Assume 
that  the  aggregate  demand  consists  of  the  sum  of  n identical  demands. 
Let  (q*  Z)  (l/n)q*.  Thus,  a fraction  of  customers  are  blacked  out 
and  they  lose  their  ability  to  consume.  Their  loss,  area  M,  is  greater 
than  L,  the  minimum  lost  consumer  surplus.  On  the  other  hand,  the 
greatest  loss  in  consumer  surplus  would  involve  area  N,  and  the  highest 
valued  uses  would  go  unserved.  A blackout  affecting  a special  set  of 
consumers  might  involve  such  a situation.  The  index  involving  lost 
quantities  does  not  necessarily  reflect  lost  values.  Engineers  use 
demand  not  served  as  a reliability  index.  Economists,  however,  would 
stress  that  the  economic  value  of  the  demand  not  served  is  the  key 
concept  for  evaluating  alternative  pricing/capacity  choices. 

As  demonstrated  by  BJ,  when  E(W)  is  maximized  with  respect  to 
price,  the  first  order  conditions  for  maximization  of  the  welfare  func- 
tion indicate  that  price  should  equal  short  run  marginal  operating  cost. 
The  optimal  price  is  lower  under  the  risk  model  than  it  is  in  the  tradi- 
tional deterministic  model  (where  p = b + p).  Given  that  the  optimal 
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Figure  2-2.  Consumer  Surplus  Versus  Demand  Not  Served 
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price  is  below  the  riskless  optimal  price,  the  law  of  demand  would  seem 
to  suggest  that  optimal  output  will  be  greater  than  equilibrium  output 
in  the  riskless  model  (Brown  and  Johnson,  1969). 

The  reason  that  prices  equal  operating  cost  can  be  clearly  ex- 
plained by  using  the  diagram  (Figure  2-3)  in  Tschirhart's  paper  (1979), 
where  four  demand  curves  are  drawn  for  four  values  of  the  disturbance 
term  (y^).  These  values  of  0.5,  0.71,  1,  and  2 enter  into  the  demand 
function  in  multiplicative  form.  Since  capacity  is  chosen  ex  ante, 
total  capacity  cost  is  essentially  a fixed  cost  (sunk  cost)  when  demand 
becomes  known.  For  demand  curve  0.5x,  price  = b is  clearly  the  welfare 
maximizing  price  since  operating  cost  b is  covered  and  consumers'  sur- 
plus is  at  a maximum.  This  price  also  applies  for  any  demand  curve  up 
to  curve  0.71x. 

For  demand  curve  2x,  if  the  lowest  valued  users  can  be  costlessly 
identified  and  not  served,  any  price  between  zero  and  p11  will  yield  the 
same  welfare.  Note  that  the  distribution  of  welfare  between  producers 
and  consumers  will  vary.  Prices  near  zero  for  demand  curve  2x  will 
require  more  rationing,  with  excess  demand  becoming  smaller  when  the 
price  is  getting  closer  to  p1'.  This  reasoning  also  applies  for  any 
demand  beyond  curve  0.71x.  Thus,  price  plays  a major  role  in  the  deter- 
mination of  service  reliability.  Tschirhart  concludes  that  in  the  short 
run,  P*  is  also  the  optimum  but  not  the  unique  price  for  all  demands 
beyond  0.71x.  Thus,  P*  was  selected  because  it  is  the  optimum  for  all 
four  cases.  The  argument  crucially  relies  on  the  assumptions  that  those 
who  value  the  service  least  are  the  first  to  be  cut  off  when  rationing 
is  necessary,  and  that  rationing  is  costless. 
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Figure  2-3.  Comparison  of  Optimal  Pricing 


19 


Meyer  (1975)  and  Crew  and  Kleindorfer  (1978)  extended  the  BJ  model 
in  order  to  rectify  this  inadequacy.  In  the  former,  a chance  constraint 
on  reliability  is  used,  while  in  the  latter,  rationing  costs  are  intro- 
duced. A chance-constrained  programming  formulation  allows  one  to 
directly  introduce  a "reliability  of  service"  or  risk  constraint  and 
analyze  the  effects  of  risk  averse  behavior  (i.e.,  demand  of  Reliable 
service)  without  the  explicit  introduction  of  a utility  function.  In 
addition,  it  provides  a formal  basis  for  analyzing  the  issue  raised  by 
BJ  (1970)  when  they  state:  "To  price  so  high  that  there  is  rarely 
(never?)  excess  demand  implies  persistent  excess  supply  (idle  capacity) 
and  an  ad  hoc  decision  by  the  public  utility  to  remove  the  risk  of 
failure.  ..."  This  approach  is  used  by  Meyer  to  obtain  prices  that 
are  sufficiently  high  to  meet  specified  standards  of  system  reliability. 
The  procedure  is  to  add  the  following  constraints  to  the  problem: 

P(xi(pi)ui  < z}  > Ei  > 0,  (2-2) 
where  P indicates  probability.  xi(p.)y.'  is  the  stochastic  demand  func- 
tion; z is  the  capacity  level.  The  value  of  the  reflects  the  strin- 
gency of  system  reliability  requirements:  a larger  e.  implies  more 
stringent  requirements.  Also  can  be  considered  as  an  LOLP  index. 
This  chance  constraint  approach  to  excess  demand  is  representative  of 
current  practice  in  some  countries.  For  instance,  planners  in  the  U.S. 
electric  power  generation  system  tend  to  use  a "one-day-in-ten  year" 
loss-of-load  probability  as  a reliability  target. 

BJ  (1970)  acknowledged  the  cost  of  excess  demand,  vis.  rationing, 
in  a reply  to  comments  by  Salkever  (1970)  and  Turvey  (1970),  but  they 
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made  no  attempt  to  add  this  cost  to  the  model . Suppose  a penalty  cost 
of  <j>  is  assessed  for  each  unit  of  excess  demand;  <f>  can  be  thought  of  as 
the  cost  of  rationing,  including  the  cost  of  ranking  consumers.  If  the 
excess  demand  only  occurs  during  the  peak  load  period,  rationing  costs 
are  a peak  load  phenomenon. 

Crew  and  Kleindorfer  (1978),  hereinafter  referred  to  as  CK,  exam- 
ined the  problems  of  optimal  pricing  and  capacity  planning  when  diverse 
plant  types  of  differing  cost  characteristics  are  available  and  when 
demand  is  stochastic.  They  concluded  that  the  existence  of  a diverse 
technology  (i.e.,  multiple  plant  types)  gives  rise  to  a system-wide 
marginal  cost  curve  which  is  upward  sloping  (this  is  a natural  conse- 
quence of  a diverse  technology),  since  plants  are  used  in  order  of  their 
equipment's  running  efficiency.  Accordingly,  marginal  cost  pricing  is 
desirable  under  peak  load  pricing  with  uncertainty.  And  on  the  other 
hand,  capacity  planning  involves  deciding  the  reliability  of  supplies. 
When  capacity  is  exceeded,  certain  'rationing  costs'  arise  in  allocating 
output  to  consumers.  The  inclusion  of  such  rationing  costs  in  the 
analysis  highlights  the  trade-off  between  losses  due  to  insufficient 
capacity  and  costs  of  additional  capacity.  Recognition  of  rationing 
costs  tends  to  further  enhance  the  desirability  of  peak  load  pricing  to 
smooth  peak  consumption  and  decrease  the  frequency  of  forced  rationing. 
That  is,  the  presence  of  rationing  costs  leads  to  an  even  more  pro- 
nounced peak  load  pricing  policy. 

CK  also  pointed  out  in  their  paper  that  the  BJ  framework  had  ne- 
glected the  important  issue  of  reliability  and  that  imposing  explicit 
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constraints  limiting  the  probability  of  excess  demand  might  signifi- 
cantly alter  their  results.  CK's  results  suggest  maximizing  the  tradi- 
tional welfare  function  subject  to  reliability  constraints.  It  is  clear 
that  the  higher  the  required  reliability  level,  the  lower  the  expected 
rationing  cost  will  be.  Thus,  the  reliability  constraint  may  be  re- 
garded as  a constraint  on  expected  rationing  costs  and  may  be  thought  of 
as  a surrogate  for  rationing  costs.  The  presence  of  such  rationing 
costs  then  induces  optimal  behavior  close  to  that  implied  by  the  usual 
deterministic  analysis.  In  particular,  peak  load  pricing  continues  to 
be  optimal  in  the  stochastic  case  since  marginal  rationing  costs  are 
higher  in  peak  periods. 

Sherman  and  Visscher  (1978),  hereinafter  referred  to  as  S V,  deve- 
loped second  best  pricing  rules  so  that  the  firm  could  remain  finan- 
cially viable.  As  Meyer  indicated  in  his  paper  (1975),  "optimal  price 
equals  marginal  (and  average  variable)  cost,  thus  a firm's  net  profits 
will  always  be  negative."  When  demand  is  stochastic,  the  expected 
welfare-maximizing  price  does  not  produce  sufficient  revenue  to  cover 
costs.  So  by  imposing  a break-even  constraint,  SV  demonstrated  how  risk 
will  influence  optimal  second  best  prices.  The  constrained  welfare 
maximizer  raises  peak  price  above  marginal  costs  by  the  greatest  amount, 
but  only  by  an  amount  sufficient  to  break  even.  When  demand  is  random, 
there  is  no  assurance  that  any  price  or  capacity  level  chosen  before 
demand  is  known  will  clear  the  market  once  demand  is  revealed.  For  if 
prices  are  set  high  enough  to  break  even  on  average,  capacity  will  go 
unutilized  whenever  demand  is  unusually  low. 

Tschirhart  (1979)  relaxed  the  implicit  assumption  that  demand  is 
independent  of  the  reliability  of  service.  The  addition  of  reliability 
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in  the  demand  function  brings  the  models  of  public  utilities  with  sto- 
chastic demand  closer  to  a position  of  being  able  to  assess  accurately 
the  welfare  implications  of  joint  pricing  and  capacity  choices.  Reli- 
ability is  an  important  aspect  of  utility  service  that  has  not  been 
fully  appreciated.  In  particular,  Tschirhart  indicated  that  an  increase 
in  price  may  improve  welfare  via  an  improvement  in  reliability,  when 
reliability  is  allowed  as  an  argument  (exogenous  variable)  in  the  demand 
function.  In  Chapter  5 we  investigate  the  cost  of  reliability  as  an 
endogenous  variable. 

Turning  from  theoretical  studies  to  more  applied  analyses  of  relia- 
bility, we  note  that  the  estimation  of  outage  cost  and  supply  cost  is 
necessary  to  achieve  optimal  decisions.  The  term  "outage  costs"  is  used 
to  encompass  all  the  economic  costs  suffered  by  society  when  the  supply 
of  electricity  is  not  perfectly  reliable.  Electric  power  supply  short- 
ages are  exhibited  in  various  ways,  ranging  from  power  surges,  frequency 
variation,  and  voltage  level  drops  (brownouts)  to  load  shedding  and 
complete  interruptions  of  supply  (blackouts).  These  effects  impose 
certain  economic  costs  on  consumers.  Mohan  Munasinghe  (1979)  developed 
a theoretical  framework  to  estimate  the  residential  outage  costs. 

The  application  of  an  optimization  technique  for  reliability  level 
requires  estimation  of  electricity  supply  costs.  Munasinghe  concludes 
that  no  generalized  production  function  model  derived  from  neoclassical 
economic  theory  could  hope  to  capture  the  detailed  characteristics  of  a 
physical  power  system.  In  particular,  sector-wide  production  functions 
based  on  cross  section  or  pooled  data  are  generally  appropriate  to 
determine  historical  trends  in  the  economics  of  electric  power  sector 
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evolution.  However,  such  functions  do  not  represent  a given  system 
accurately,  because  the  estimates  remain  oriented  in  the  past  rather 
than  the  present  or  future  (i.e.,  they  do  not  incorporate  current  input 
prices  or  units  coming  on-line).  An  engineering  model  used  in  the 
design  and  planning  of  a system  incorporates  future  system  supply  costs 
when  the  reliability  level  is  to  be  optimized.  Of  course  the  process  of 
system  planning  and  determining  electricity  supply  costs  is  usually 
complicated  by  several  practical  difficulties,  including  resolving 
uncertainties  in  the  demand  projections,  errors  in  the  cost  estimates, 
technological  changes,  construction  delays,  and  environmental 
constraints. 

2.3  Empirical  Studies  of  Reliability 

The  optimum  level  of  reliability  is  easier  to  define  theoretically 
than  to  measure  in  practice.  Its  computation  requires  two  sets  of  data: 

1.  A set  which  shows  how  the  incremental  costs  (3C)  respond  to 
additional  reliability  of  supply  (3R),  that  is, 

2.  Another  set  of  data  showing  to  what  extent  consumer  loss  (9L) 
is  spared  and  alleviated  through  increased  reliability  (3R), 
that  is,  - (this  value  is  negative  because  increased  relia- 
bility causes  less  loss) . 

No  previous  empirical  studies  have  established  these  components 
with  accuracy.  Appropriate  data  are  unavailable  at  either  the  micro 
(individual  firm)  or  sectoral  level.  Tschirhart  (1980)  indicated  that 
empirical  work  on  demand  estimation  with  reliability  as  an  argument 
would  be  difficult  in  most  countries  since,  historically,  reliability 
has  been  on  the  order  of  99  percent  with  little  variation.  His  analysis 
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pointed  out,  however,  that  there  may  be  welfare  gains  available  through 
decreases  in  reliability  because  over-reliable  service  is  simply  too 
costly.  Consumers,  if  given  a chance,  may  opt  for  lower  reliabilities 
if  accompanied  by  lower  prices.  A reliability  of  75  percent  would 
undoubtedly  be  unacceptably  low,  but  a "one-day-in-one-year"  (99.73 
percent)  or  a "one-day-in-f ive-year"  (99.95  percent)  loss  of  load  proba- 
bility may  be  a better  target  than  the  "one-day-in-ten-year"  (99.97 
percent)  estimate. 

Sant  (1979)  suggested  that  cost  functions  be  modeled  to  include 
both  yearly  output  and  the  temporal  distribution  of  that  output  (i.e., 
load  factor).  Otherwise,  he  suggested,  the  estimate  of  electricity 
costs  would  be  biased  due  to  generation  mix  effects.  He  still  failed  to 
stress  that  the  reliability  of  electricity  should  also  be  included  in 
the  cost  equation,  although  the  load  factor  could  be  viewed  as  a proxy 
for  how  reliability  would  be  affected  by  peak  KW  demands. 

Michael  Telson  (1975)  examined  the  cost  of  providing  an  incremental 
amount  of  production  through  installing  an  additional  amount  of  genera- 
tion to  a system  which  was  built  at  the  loss  of  load  probability  cri- 
terion of  one-day-in-ten-years,  and  found  that  society  would  be  paying 
at  least  208  times  the  revenue  of  the  additional  energy  produced  by  the 
incremental  megawatt.  He  also  examined  the  welfare  losses  caused  by 
not  having  that  energy  and  concluded  that  we  could  at  most  attribute  a 
factor  of  forty  to  this  value.  Thus,  his  calculations  suggest  that 
one-day-in-ten-year  loss  of  load  probability  criterion  is  too  conserva- 
tive. He  recommends  reducing  generation  system  reliability  target 


25 


levels  to  perhaps  a f ive-day-in-ten-year  criterion.  It  is  also  inter- 
esting to  note  that  the  United  States  seems  to  have  had  the  highest 
electric  service  reliability  experience  in  the  world.  He  also  quoted 
the  U.S.  experience  as  99.98  percent  reliability  in  terms  of  energy 
sales  and  a 1.75-hour  loss  per  customer  per  year.  The  equivalent  fig- 
ures for  the  next  highest  countries,  France  and  Great  Britain,  were 
99.93  percent  and  six  hours  and  99.80  percent  and  17.5  hours,  respec- 
tively. As  the  cost  of  reliability  rises,  it  will  become  more  attrac- 
tive to  reduce  reliability  target  levels  and  to  devise  pricing  and/or 
rationing  schemes  to  reduce  the  peak  loads  which  account  for  a great 
part  of  an  electric  system's  unreliability,  unless  the  benefits  of 
reliability  increase  by  the  same  amount.  However,  note  that  modern 
microprocessor  metering  technology  and  interruptible  rate  structures 
make  it  possible  to  reduce  the  capacity  requirement  for  given  stochastic 
demands . 

Crew  and  Kleindorfer  (1979)  used  ordinary  least-squares  (OLSQ)  to 
estimate  a reliability  equation.  The  proposed  equation  form  was  as 
follows  (all  variables  are  measured  in  the  current  period  t except  for  s 
which  was  lagged  k periods): 

Rt  " f(St-k’  Lt’  It’  Et’  Nt’  Ht’  Ct’  Gt’  V‘ 

Data  from  39  utility  companies  covered  a period  from  1971  to  1978,  and 

consisted  of  the  following: 

S = Reported  rate  of  return  on  total  assets 
L = Load  factor  (a  measure  of  utilization  of  a system) 

I = Imports  of  power 
E = Exports  of  power 
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N = Proportion  of  nuclear  capacity 
H = Proportion  of  hydro  capacity 
C = Proportion  of  fossil  capacity 
G = Proportion  of  internal  combustion  capacity 

A = Average  plant  size 

r = Capacity  - peak 
Capacity 

This  is  the  only  published  paper  to  discuss  the  estimation  of  the  relia- 
bility equation.  Also,  only  one  reliability  index,  reserve  margin,  was 
used  in  this  analysis. 

Their  concern  was  mainly  with  the  effects  of  regulation  on  poten- 
tial overcapitalization  and  reliability.  Crew  and  Kleindorfer  concluded 
that  the  inclusion  of  reliability  in  the  model  offers  no  straightforward 
way  to  test  the  AJ  effect  under  uncertainty.  Their  use  of  ex  post 
reserve  margin  as  the  reliability  index  severely  limits  their  analysis. 
Furthermore,  by  definition,  installed  capacity  is  used  to  calculate 
reserve  margin;  however,  installed  capacity  is  nonstochastic.  Both 

demand  and  supply  are  stochastic  variables  in  the  calculation  of  LOLP 
and  e(DNS).  Taking  another  approach  to  the  problem  of  cost,  James 
Griffen  (1977)  modeled  a long  run  electric  power  generation  cost  func- 
tion with  pseudo  data.  He  used  a process  model  to  generate  the  data. 
The  electric  power  process  model  used  in  his  paper  consisted  of  73 
linear  equations  with  98  process  activities.  The  process  model  is  a 
linear  programming  problem,  with  the  following  translog  cost  function 
model  as  follows: 

, 8 8 8 

In  C - lna0  + In  Q + ^ o^.p.  + ^ • E 8ij (In  p^ln  p. , 
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where  C is  total  cost,  Q is  output,  and  p^.  refers  to  the  eight  input 
prices.  Surprisingly , he  did  not  include  any  reliability  variables  in 
the  cost  model . 

According  to  Griffen,  his  pseudo  data  approach  offered  numerous 
advantages  compared  to  conventional  time  series  particularly  in  that 
they  avoid  mul ticol 1 ineari ty , a limited  sample  range,  and  inadequate 
technical  and  environmental  detail.  He  noted  that  the  pseudo  data 
approach  is  subject  to  certain  caveats  and  limitations.  The  value  of 
the  pseudo  data  approach  ultimately  rests  on  the  quality  and  complete- 
ness of  the  underlying  process  model,  which  is  typically  difficult  for 
non-engineers  to  validate.  Pseudo  data  are  also  subject  to  both  mea- 
surement errors  and  errors  regarding  assumed  behavioral  parameters. 

Smith  and  Vaughan  (1979)  pointed  out  that  the  results  derived  from 
Griffen's  exercise  are  only  as  good  as  the  information  in  the  underlying 
process  analysis  model,  including  the  detail  with  which  the  technology  is 
represented.  They  used  three  process  analysis  models  to  analyze  iron 
and  steel  production.  Their  findings  suggest  that  Griffen's  warning 
concerning  the  possible  limitations  of  pseudo  data  was  well  founded. 

According  to  Maddala  and  Roberts  (1981),  it  is  not  clear  that  the 
way  Griffen  went  about  condensing  the  technical  information  in  the 
process  model  is  very  useful  because  Griffen's  pseudo  data  approach 
deals  with  an  attempt  at  a single-equation  summarization  of  data  gener- 
ated from  a complex  process  analysis  model  (i.e.,  a linear  programming 
problem  was  set  up  and  solved).  However,  they  did  agree  to  Griffen's 
claim  that  one  cannot  study  the  effects  of  new  technology  and  environ- 
mental policies  with  the  help  of  time-series  data,  and  that  one  has  to 
look  at  more  micro-engineering  data  to  get  an  idea  of  these  effects. 
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In  this  study,  due  to  lack  of  significant  and  sufficient  empirical 
data,  we  will  simulate  a process  model  to  generate  pseudo  data  for 
empirical  tests.  This  approach  promises  to  yield  insights  into  system 
reliability.  Before  discussing  the  particular  prices,  technological 
demand  parameters,  and  the  respective  stochastic  components,  we  need  to 
have  a clear  idea  of  how  various  reliability  indices  are  related.  Such 
detailed  discussion  will  be  presented  in  Chapter  3. 

2.4  Directions  for  Research 

One  of  the  management  options  which  we  believed  needed  further 
research  was  the  use  of  pricing  policies.  Pricing  was  considered  to  be 
an  important  facet  of  shortage  prevention  and  management.  Research 
should  cover  the  effectiveness  of  various  pricing  schemes  and  potential 
draw-backs  in  regulatory-determined  rate  structures  for  utilities.  The 
equipment  and  administrative  costs  of  implementing  and  managing  alter- 
native price  schemes  were  also  considered  an  area  in  need  of  further 
research.  For  example,  new  metering  devices  would  be  required  for 
residential  customers  for  the  introduction  of  time-of-day  pricing.  The 
costs  of  such  meters  need  to  be  taken  into  account  if  such  rate  struc- 
tures are  to  be  used  in  a cost-effective  manner. 

In  addition,  models  should  be  developed  to  assess  the  socio- 
economic impacts  of  an  electricity  shortage.  The  development  of  appro- 
priate criteria  for  optimizing  reliability  is  also  an  important  issue  in 
need  of  more  research.  Finally,  sound  process  models  are  needed  to 
generate  pseudo  data  for  facilitating  current  research  and  to  provide  a 
guideline  for  future  data  collection. 


CHAPTER  3 

RELIABILITY  INDICES  AND  COST  CHARACTERISTICS 
FOR  ELECTRIC  UTILITIES 

The  term  reliability  is  used  in  an  engineering  sense  to  indicate 
quality  of  constancy  of  service.  A number  of  indices  have  been  intro- 
duced in  reliability  theory  to  characterize  reliability,  such  as  Effec- 
tive Reserve  Margin  (RM),  Loss  of  Load  Probability  (LOLP),  Expected 
Demand  not  Served  (e( DNS) ) , Expected  Energy  not  Served  (e(ENS) ) , Ex- 
pected Loss  of  Load  (XLOL),  Frequency  and  Duration  (FAD),  and  Safety 
Margin  (S).  Just  as  a four-firm  concentration  index  cannot  fully  repre- 
sent all  the  dimensions  of  market  power,  no  single  reliability  index  is 
likely  to  be  an  ideal  measure  of  the  multi-dimensional  quality  variable 
called  "reliable  service."  For  instance,  we  find  electricity  systems 
planners  moving  from  the  traditional  reserve  margin  index  to  a loss-of- 
load-probability  index  for  investment  planning  purposes.  Given  the 
amount  of  resources  involved,  the  topic  is  important  both  to  utility 
companies  and  regulators. 

This  chapter  presents  mathematical  definitions  and  concise  analyses 
of  various  reliability  indices,  including  the  advantages  and  disadvan- 
tages of  each  index.  Their  links  to  cost  are  especially  important  since 
when  evaluating  the  reliability  of  a number  of  alternative  power  plans, 
it  is  necessary  to  estimate  the  cost  associated  with  a particular  plan. 
Thus,  we  will  also  review  the  cost  characteristics  of  electric  utilities 
in  the  later  sections  of  this  chapter. 
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3.1  Reliability  Indices 

3.1.1  Effective  Reserve  Margin  (RM) 

Perhaps  the  most  commonly  used  index  in  the  press  is  the  reserve 
margin,  although  reported  numbers  are  ex  post  '(after  stochastic  demand 
and  supply  have  been  realized).  This  index  in  its  ex  ante  form  can  be 
defined  as 

IC  - Le 


RM  = 


IC 


(3-1) 


where 


IC  = installed  capacity,  and 

Le  = effective  load, 

= L + Z L . . 

01 

Again,  L = expected  load  (ex  ante),  and 


Loi=  the  random  outage  load  of  the  ith  unit. 

The  rationale  for  this  calculation  is  that  the  firm  should  provide  for 
such  factors  as  unit  forced  outages,  scheduled  outages  for  maintenance, 
and  load  forecast  uncertainty,  when  relating  the  effective  load  to  the 
available  capacity.  Sometimes  the  ex  post,  or  realized  reserve  margin, 
is  used  as  an  index  of  reliability,  but  this  version  of  the  index  is 
invalid  as  a planning  tool. 

3.1.2  Loss  of  Load  Probability  (LOLP) 

LOLP  is  defined  as  the  probability  that  effective  load  exceeds 
installed  capacity  (IC)  of  a particular  proposed  generation  system.  At 
present,  the  LOLP  criterion  is  commonly  used  for  generation  planning. 
Since  system  failure  is  most  often  defined  in  terms  of  inability  to  meet 
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the  daily  peak  load,  the  LOLP  is  usually  computed  not  as  a strict  proba- 
bility but  rather  as  an  expectation:  This  index  may  be  interpreted  as 

the  average  number  of  days  over  a long  period  during  which  the  daily 
peak  demand  is  expected  to  exceed  the  available  generating  capacity. 
For  example,  the  probability  that  the  load  will  exceed  the  installed 
capacity  of  100  MW  of  the  generation  system  is  0.02981;  that  is,  LOLP  = 
0.02981.  Because  the  number  0.02981  is  unclear,  it  is  common  to  convert 
the  LOLP  into  so  many  hours  per  week  or  days  per  year,  because  the  LOLP 
can  be  interpreted  as  the  percentage  of  time  the  load  is  expected  to 
exceed  the  installed  capacity.  Here,  the  load  is  interpreted  as  in- 
cluding outages  or  down-time  for  the  capacity. 

LOLP  (in  hours/week)  = 0.02981  x 168 

= 5.008  hours/week, 
or 

LOLP  (in  days/year)  = 0.02981  x 365  days/year 

= 10.88  days/year. 

The  target  LOLP  frequently  used  for  planning  the  United  States  is  one- 
day-in-ten-years,  but  the  rationale  for  this  target  is  unclear.  As  was 
noted  earlier,  Telson,  among  others,  is  critical  of  this  goal  (1975). 

3.1.3  Expected  Demand  Not  Served  ( e( DNS) ) 

In  order  to  understand  the  next  index,  expected  demand  not  served, 
we  need  to  introduce  the  effective  load  probability  density  function. 
Suppose  a continuous  random  variable,  effective  load  (denoted  by  X),  has 
a normal  distribution,  with  a mean  of  50  MW  and  a standard  deviation  of 
15  MW.  Assume  that  we  have  70  MW  of  installed  capacity.  What  is  the 
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probability  that  the  effective  load  X will  be  more  than  70  MW  (installed 
capacity)?  By  definition,  this  probability  is  LOLP.  Also,  this  proba- 
bility is  the  area  under  the  normal  probability  distribution  more  than 
70  MW,  as  shown  in  the  shaded  area  A of  Figure  3-1 (A) . 

The  first  step  in  finding  the  area  A is  to  calculate  the  Z-score 

corresponding  to  the  measurement  70  MW.  We  calculate 

Z = X - u 
a 


Thus,  the  measurement,  70  MW,  is  1.33  standard  deviations  away  from 
the  mean,  50  MW.  By  looking  at  the  Z-table,  we  find 
p( 50  < x < 70)  = 0.4082 
= area  B. 

The  desired  probability  is 
p(Z  > 1.33)  = 0.5  - B 

= 0.5  - 0.4082 
= 0.0918  = LOLP  at  70  MW. 

We  can  repeat  the  above  procedure  to  derive  the  effective  load  proba- 
bility distribution  as  shown  in  Figure  3- 1 ( B ) . For  example,  the  proba- 
bility that  the  effective  load  will  exceed  50  MW  is  0.5. 

Expected  demand  not  served  (e(DNS))  is  beginning  to  see  widespread 
use  because  of  its  obvious  economical  significance,  and  is  defined  as 
the  expected  value  of  instantaneous  KW  demand  in  excess  of  installed 
capacity,  IC.  In  order  to  explain  this  concept,  we  refer  to  Figure  3-2. 
By  definition,  if  (e(DNS))  = 0,  all  units  are  100  percent  reliable  and 
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f (x ) A.  Effective  load  with  a normal  distribution 


Figure  3-1.  Random  Loads  and  the  Effective 
Load  Probability  Distribution 
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Figure  3-2.  Typical  Effective  Load  Probability 
Di stribution 
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capacity  is  greater  than  the  maximum  expected  MW  demand.  The  mathemati- 
cal definition  of  (e(DNS))  is 

o(DNS)  - F(Le)dLe,  , (3-2) 

which  is  the  area  under  the  curve  in  the  above  figure  from  Le  = IC  to 
0 + IC,  as  indicated. 

3.1.4  Expected  Energy  Not  Served  (e(ENS)) 

If  Equation  3-2  is  multiplied  by  T,  the  time  period  represented  by 
F(Le),  the  (e(ENS))  results: 

(e(ENS) ) = T / l+elCK  F(Le)  dLe  (3-3) 

This  concept  is  more  closely  related  to  economic  rather  than  engineering 
features  in  that  the  amount  (although  not  the  dollar  value)  of  consump- 
tion is  being  measured.  The  unit  of  measurement  for  e(DNS)  is  MW,  while 
the  measure  for  e(ENS)  is  MWH. 

3.1.5  Expected  Loss  of  Load  (XLOL) 

The  LOLP  is  often  supplemented  by  the  expected  loss  of  load  index, 
which  indicates  the  expected  magnitude  of  the  unsupplied  load,  given 
that  a failure  has  occurred;  therefore,  XLOL  is  a conditional  expecta- 
tion. That  is, 

[load  - supply/load  - supply  > 0]  = e(DNS)/L0LP  . (3-4) 

3.1.6  Frequency  and  Duration  (FAD) 

The  expected  frequency  or  mean  recurrence  time  between  outages  and 
their  expected  duration  represent  two  other  reliability  indices.  These 
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indices  are  used  widely  in  transmission  and  distribution  system  plan- 
ning, but  rather  infrequently  in  generation  studies,  and  can  be  cate- 
gorized as  follows: 

i)  Frequencies,  such  as  the  average  number  of  failures  per  unit 
time,  and 

ii)  Mean  duration,  such  as  the  mean  time  to  the  first  failure  and 
the  mean  time  between  failures. 

3.2  Comparisons  of  Indices 

Table  3-1  summarizes  commonly  used  reliability  indices.  The  defi- 
nitions and  units  of  measurement  serve  to  illustrate  the  wide  range  of 
concepts  used  by  practitioners.  Several  general  points  about  measures 
of  reliability  should  be  noted.  First  of  all,  each  of  the  indices  has 
its  strengths  and  weaknesses,  and  cannot  individually  provide  a complete 
description  of  outage  patterns.  For  example,  although  LOLP  is  a measure 
of  the  probability  of  system  failure,  the  LOLP  does  not  indicate  the 
expected  magnitude  of  the  load  lost  during  such  a contingency.  In 
contrast,  the  XLOL  provides  the  expected  magnitude,  but  neither  the 
probability  nor  the  average  duration  of  an  outage.  Like  the  LOLP,  the 
FAD  measure  also  does  not  indicate  the  mean  size  of  the  failure,  but 
because  the  FAD  indices  are  specified  at  the  load  point  or  particular 
point  of  consumption,  the  number  of  affected  customers  may  be  deduced. 

Second,  the  methods  of  computation  and  their  interpretation  may 
vary.  Although  they  are  often  related  in  well-defined  ways,  the  indices 
are  generally  neither  comparable,  nor  consistent:  they  could  yield 
different  rankings  of  systems  by  reliability  level.  Similarly,  rankings 
by  cost  over  the  planning  horizon  can  be  expected  to  differ,  although 
establishing  a baseline  for  making  comparisons  is  no  simple  task. 
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Third,  the  estimated  indices  may  be  much  more  significant  in  a 
relative  rather  than  an  absolute  sense.  Reliability  indices  are  usually 
computed  as  predictors  of  long  run  average  system  performance,  but  they 
do  not  reflect  actual  experience  well.  That  is,  the  difference  between 
actual  and  predicted  values  can  be  substantial.  The  range  depends  on 
the  extent  of  potential  forecast  errors,  in  terms  of  unit  failures  and 
electricity  demand. 

Fourth,  the  LOLP  index  is  more  significant  in  the  mind  of  the 
system  planner,  but  from  the  consumer's  standpoint,  the  FAD  is  more 
significant.  The  FAD  index  is  easier  to  understand,  and  production  and 
consumption  activities  are  directly  related  to  the  duration  and  fre- 
quency of  outage.  Given  particular  FAD  patterns,  the  e(DNS)  index  may 
be  of  greatest  interest  to  the  economist,  since  it  raises  two  important 
economic  issues:  the  opportunity  cost  of  outage  and  the  supply  cost  of 

eliminating  outage. 

Lastly,  the  implications  of  alternative  indices  warrant  attention. 
Let  us  begin  by  reproducing  the  right  corner  of  Figure  3-3  (see  Figure 
3-4).  Again,  the  shaded  area  represents  the  e(DNS),  and  the  intercept 
of  the  vertical  axis  represents  the  corresponding  LOLP.  We  can  consider 
three  illustrative  cases. 

Case  A.  Both  Curve  A and  Curve  B in  Figure  3-4  represent  different 
magnitudes  of  e(DNS).  That  is,  e(DNS)  of  Curve  A is  greater  than  the 
e(DNS)  of  Curve  B.  However,  Curve  A and  Curve  B have  the  same  LOLP 
denoted  by  point  L.  Can  we  claim  that  two  curves  have  the  same  relia- 
bility? In  order  to  answer  this  question,  it  is  necessary  to  introduce 
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Figure  3-3.  LOLP  and  e(DNS) 
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Figure  3-4.  Same  LOLP  vs.  Different  e(DNS) 
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a new  concept,  "generalized"  reliability,  denoted  by  R*.  R*  is  a func- 
tion of  all  reliability  indices.  That  is  R*  (Rj,  R£,  . . . , R ).  For 
example  R^  = LOLP,  R2  = e(DNS),  etc. 

For  the  convenience  of  comparison,  we  define 
★ 

R^  = Generalized  reliability  preference  of  Curve  A,  and 
★ 

Rg  = Generalized  reliability  preference  of  Curve  B. 

In  this  particular  example 
★ * 

Rg  > R^:  B is  more  reliable  than  A,  where 

R^l=  LOLP  of  Curve  A, 

Rgl=  LOLP  of  Curve  B, 

R^2=  e(DNS)  of  Curve  A,  and 
Rg2=  e(DNS)  of  Curve  B. 

"k  k 

In  Figure  3-4,  RA(RAj,  ^A2^  < ^B^Bl’  RB2^‘  This  1S  simP^  because  R^ 
= Rgp  Rg2  < R^2*  Thus,  the  answer  for  Case  A is  that  Curve  B is  more 
reliable  than  Curve  A. 

Case  B . Both  Curve  A and  Curve  B,  shown  in  Figure  3-5,  represent  the 
same  magnitude  of  e(DNS).  However,  they  have  different  LOLP  indices, 
denoted  by  point  L and  point  M respectively.  Again,  can  we  claim  that 
both  Curve  A and  Curve  B have  the  same  reliability;  by  the  same  decision 

★ k 

criterion:  Rftl  > Rgl  and  R^  = Rg2»  so  Rft(RA1,  RA2)  < Rg(RB1>  RB2'*  In 
other  words,  Curve  B is  more  reliable  than  Curve  A,  if  these  are  the 
only  two  dimensions  of  reliability  that  matter.  If  we  bring  XLOL  into 
the  above  function  (and  drop  LOLP),  we  have 
XLOLg  > XL0La, 

since  XLOL  = eDNS/LOLP.  Thus,  XL0LA  is  the  area  under  the  A curve 
divided  by  the  LOLP:  Since  the  areas  are  the  same  but  LOLP^  is  greater 
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Figure' 3-5.  Different  LOLP  vs.  Same  e(DNS) 
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than  LOLPg,  XLOL^  is  less  than  XLOLg.  Thus,  we  would  conclude  that 
curve  A is  more  reliable  than  curve  B. 

Case  C.  Even  though  both  Curve  A and  Curve  B,  shown  in  Figure  3-6,  have 
the  same  magnitude  of  LOLP  and  e(DNS),  it  is  still  inappropriate  to 
claim  that  both  Curve  A and  Curve  B have  the  same  reliability.  For 
instance,  e(DNS)  of  curve  B might  be  caused  by  just  a one-shot  loss. 
However,  e(DNS)  of  curve  A is  accumulated  by  many  small  losses.  Thus, 
both  curves  probably  represent  different  outage  frequency  and  duration. 
Also,  both  curves  might  have  different  cost  implications.  For  example, 
the  power  generating  system  of  curve  B could  consist  of  two  generators 
(80MW  and  20MW),  and  the  power  generating  system  of  curve  A could  con- 
sist of  five  generators  (5MW,  10MW,  15MW,  25MW,  and  45MW).  Since  the 
generating  systems  have  different  capacity  mixes,  they  have  different 
cost  structures,  even  if  they  have  the  same  total  capacity.  The  latter 
might  cost  a little  more  due  to  more  generators  in  the  capacity  mix, 
resulting  in  a reduction  in  scale  economies. 

In  conclusion,  each  of  the  reliability  indices  has  its  own 
strengths  and  weaknesses  and  cannot  individually  give  a complete  de- 
scription of  power  system  reliability.  They  could  yield  different 
rankings  of  the  system  by  reliability  level.  We  need  to  consider  all 
indices  when  evaluating  power  system  plans,  and  the  weight  assigned  to 
each  index  in  the  generalized  reliability  function  should  also  be  taken 
into  consideration.  From  the  standpoint  of  public  policy,  the  weights 
ought  to  reflect  the  implications  of  the  indices  for  the  economic  value 
of  lost  consumption.  Chao  (1983)  takes  a step  in  this  direction  using 
some  specific  consumer  surplus  loss  functions. 


LOLP 


MW 


Figure  3-6.  Same  LOLP  vs.  Same  e(DNS) 


45 


3.3  A Stylized  Model  of  Utility  Costs 

Capacity  costs,  energy  costs,  and  timing  of  additions  are  three 
major  elements  which  influence  the  overall  cost  of  expansion  plans.  For 
simplicity,  we  will  use  a highly  stylized  model  of  an  electric  utility 
to  generate  a set  of  data  reflecting  the  engineering/economic  trade-off 
relating  short-run  production  costs  and  reliability.  The  numbers  are 
meant  to  illustrate  production  constraints  and  to  provide  a data  base 
for  testing  the  applicability  of  traditional  econometric  cost  functions. 
The  total  annualized  capacity  cost  of  unit  3-j  is  defined  as  follows: 

3j  = Ficciki,  and  (3-5) 

= 3*(1  >000)  = (F.cc.k.Jd.OOO)  ^ (3-6) 

1 Q Q ’ 

P 

where 

cc.j  = capacity  acquisition  cost  per  MW  of  unit  i, 

k.j  = capacity  in  MW  of  unit  i, 

Qp  = peak  demand  (KW  x peak  hours), 

cc^  ‘ k.j  = total  actual  acquisition  cost  ($/MW), 

F.j  = fixed  charge  rate  (covers  taxes,  depreciation,  and  cost 
of  capital), 

3*  = total  annualized  capacity  cost  ($/MW),  and 

3* 

3 = -fr  = capacity  cost  per  peak  KWH. 
qP 

By  definition,  the  fixed  charge  rate  associated  with  investment  in  unit 
i is  the  annual  expense  associated  with  the  investment  expressed  as  a 
percentage  of  initial  investment  (l'.e.,  cc.k.).  It  reflects  taxes, 
depreciation,  and  the  allowed  return  on  investment. 

The  energy  cost  associated  with  unit  i can  be  defined  as 
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where 

f.j  = fuel  cost  per  KWH, 

= operation  and  maintenance  cost  per  KWH,  and 
bi  = energy  cost  per  KWH  of  unit  i. 

Because  future  fuel  costs  are  very  difficult  to  predict,  utilities 
must  consider  a range  of  fuel  costs  and  fuel  types  in  evaluating  expan- 
sion alternatives.  Typical  operating  and  maintenance  costs  for  various 
plant  types  are  obtained  through  engineering  studies. 

We  can  easily  define  total  production  cost  in  terms  of  equations 
3-6  and  3-7: 

TC  = 3Qp  + b(Qo  + Qp) , (3-7) 

where 

Qp  = peak  demand,  and 
Q0  = off-peak  demand. 

Both  sides  of  Equation  3-7  can  be  divided  by  the  expected  value  of  the 
annual  energy  produced  by  unit  i,  denoted  by  e(E.),  then  the  total  cost 
might  be  expressed  in  mil s/KWH  or  $/MWH,  which  is  often  convenient  when 
comparing  different  plants  types.  Engineers  define  e(E.)  as  follows: 
e(E.)  = 8,760  e(CF.)k 
- 8,760  CF . k , 

where 

CF.  = capacity  factor. 

Clearly,  the  timing  of  investment  via  the  addition  of  a new  unit  can 
drastically  change  the  operation  of  existing  units,  thus  changing  total 


system  cost. 
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3.3.1  Cost  (capacity  cost  and  fuel  cost)  and  Mix  of 
Generating  Plant 

For  simplicity,  only  two  kinds  of  generating  units  are  considered: 
gas  turbines  and  base  load  steam  plants.  Figure  3-7  is  a graph  of  the 
(annualized)  total  cost  per  megawatt  of  installed  capacity,  that  is, 
capital  cost  plus  operating  cost,  plotted  as  a function  of  the  numbers 
of  hours  of  operation.  Gas  turbines  have  lower  capital  costs  (repre- 
sented by  the  intercept  3^),  but  higher  fuel  costs  (as  indicated  by  the 
slope,  bj,  of  the  operating  cost  curve)  than  the  steam  plant.  The  cost 
characteristic  of  these,  units  has  been  simplified  for  convenience  of 
analysis  by  assuming  fixed  values  of  $£,  b^,  and  As  the  upper 
diagram  shows,  from  the  economic  point  of  view,  gas  turbines  are  more 
expensive  than  base  load  units  if  they  are  to  be  used  more  than  H hours 
a year.  The  lower  diagram  for  load  curve  A indicates  that  it  would  be 
economical  to  serve  the  bottom  OF  megawatts  of  load  (i.e.,  the  base 
load)  with  steam  units  and  the  upper  FG  megawatts  of  load  (i.e.,  the 
peak  load)  with  gas  turbines.  For  load  curve  B,  the  choice  rule  for 
capacity  mix  is  changed. 

In  practice,  there  are  many  kinds  of  generating  units  including 
nuclear,  coal  and  oil-fired  steam,  gas  turbines,  and  conventional  and 
pumped-storage  hydro.  The  units  have  characteristics  varying  by  type 
and  by  size.  All  would  have  to  be  taken  into  account,  when  generation 
additions  were  being  considered  to  meet  future  load  growth.  However, 
the  basic  principle  of  selecting  units  with  lower  capital  costs  but 
higher  operating  costs  for  peaking  purposes  generally  holds  true. 
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Figure  3-7.  Cost  Relationships 
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3.3.2  Average  Costs  of  Base  Load  and  Peak  Load  Energy 
It  must  not  be  thought  that  every  energy  kilowatt-hour  costs  the 
same  to  produce,  even  if  generated  by  the  same  plant.  This  can  be  seen 
from  the  following  purely  illustrative  example.  Assume  that  the  load  is 
supplied  by  a single  plant  of  750  MW  capacity.  Further  assume  that  the 
generation  costs  from  this  plant  are  $10.5  million  per  annum  (capacity 
cost)  plus  2.5  mills  per  KWH  (variable).  The  load  duration  curve  may  be 
conceived  of  as  consisting  of  a very  large  number  of  thin  horizontal 
elements  such  as  that  shown  (not  to  scale)  in  Figure  3-8.  As  the  thick- 
ness of  this  element  is  one  kilowatt,  it  bears  1/750, 000th  of  the  capa- 
city cost  ($10.5  million)  or  $14/KW  per  annum.  The  energy  generated  by 
this  particular  kilowatt  during  the  year  will  be  h kilowatt-hours,  where 
h is  the  breadth  of  the  element  in  hours.  Hence  the  capacity  costs  to 
the  element  may  be  expressed  as 

Capacity  cost  per  KWH  in  mills  = 14  -x  ^ ?-QQ0-  , ($1  = 1,000  mills). 

To  this  must  be  added  the  variable  cost  of  2.5  mills  per  KWH  so  that  the 
average  costs  of  generating  this  element  of  energy  will  be 
( 14^-°--  + 2.5)  mills  per  KWH. 


Assume  that  the  particular  element  illustrated  in  Figure  3-8  corresponds 
to  a load  of  440  MW  at  which  level  the  value  of  h is  2,800  hours.  The 
average  cost  of  generating  this  element  of  energy  will  therefore  be 


AFC  + AVC  = AC 
14,000 


2,800 


+ 2.5  = 7.5 


mi  1 s/KWH, 


where 

AFC  = average  fixed  cost,  and 
AVC  = average  variable  cost. 
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Figure  3-8.  Annual  Duration  Curve  Figure  3-9.  Average  Cost  for  Various  Loads 
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This  figure  is  represented  by  p in  Figure  3-9,  which  shows  on  the  hori- 
zontal axis  the  average  costs  of  each  kilowatt  element  of  load  from  the 
base  load  to  the  extreme  peak.  Note  that  the  value  of  h is  getting 
smaller  as  the  load  moves  from  base  load  to  the  peak  load. 

In  practice,  the  load  will  be  supplied  by  several  plants.  From  the 
standpoint  of  efficiency,  it  is  desirable  that  the  base  load,  the  inter- 
mediate loads,  and  the  peak  loads  be  supplied  by  different  types  of 
plants  with  different  cost  characteristics. 


3.3.3  Numerical  Calculations  of  Capacity  Cost  (3) 

In  Figure  3-10,  we  assume  that  the  total  actual  acquisition  cost 
(cc.jk.j)  is  $109. 5/MW  and  the  fixed  charge  rate  (F^)  is  10  percent. 

By  Equation  (3-5) 


3+  = 10%  x $109.5  = $10. 95/MW. 

Further  assume  under  uniform  pricing  that  peak  KWH  demand  per  day  (12 
hrs.  peak,  12  hrs.  off-peak)  is  1,000  KWH,  and  off-peak  KWH  demand  per 
day  = 500  KWH.  Thus,  1,000  KWH/12  hrs.  is  needed  for  capacity  (i.e., 

83  1/3  KW  to  meet  hourly  peak  demand  per  day).  Furthermore,  annual  KWH 
peak  consumption  (Q  ) = 83  1/3  KW  x 12  hrs.  x 365  days  = 365,000  KWH, 
and  annual  KWH  off-peak  consumption  (Q  ) = 41  2/3  KW  x 12  hrs.  x 365 
days  = 182,500  KWH.  In  other  words,  annual  total  consumption  is  547,500 
KWH/yr. 

Capacity  cost  per  peak  KWH  can  be  derived  from  Equation  (3-6), 

Bi  i 365.000  = 3,t/KWH  (peak)' 


Also,  average  capacity  charge  per  KWH  under  uniform  pricing  is 


- - j 
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Figure  3-10.  Uniform  Pricing  vs.  Peak  Load  Pricing 
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$10,950/547,500  ( i . e . , 2<£/KWH).  Note  that  the  shaded  area  in  Figure 
3-11  shows  short  run  misallocations  (efficiency  losses)  under  uniform 
pricing: 

(i)  TXB  represents  the  loss  due  to  under-pricing  for  peak 
customers,  and 

it 

(ii)  AAKS  represents  the  loss  due  to  over-pricing  for  off-peak 
customers . 

3.3.4  Reliability  and  Peak  Load  Pricing 

As  shown  in  Figure  3-10,  peak  quantity  demanded  (1,000  KW  per  12 
hours)  under  uniform  pricing  has  been  reduced.  The  reduction  of  peak 
demand  implies  the  increase  of  reliability  and  the  decrease  of  outage 
rate.  Also,  there  are  two  welfare  gains.  One  is  capacity  saving 
(greater  than  consumer  valuation  of  the  output),  represented  by  ATXB. 
The  other  is  the  increase  of  consumer  surplus  (which  is  greater  than  the 
increase  in  fuel  costs),  represented  by  AAKS. 

3.3.5  Plant  Mix  under  Diverse  Technology 

Regarding  the  optimal  prices.  Crew  (1979)  reveals  the  result  that 
prices  are  set  at  marginal  cost.  (See  notes  for  detailed  mathematical 
proofs.)  There  are  two  major  conclusions  regarding  optimal  prices  and 
optimal  plant  mix.  Under  diverse  technology 

(i)  optimal  price  (p^  = off  peak  price,  p^  ~ Pea^  price) 

?2  ~ ^2  + ^2 5 

Pj  = (2b1  + e:)  - (b2  + 32), 

where 

bi  < Pi  < b2  < b2  + B2  = P2  < bi  + $i  • 

The  introduction  of  a more  diverse  technology  leads  to  lower 
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1 


1. 


Figure  3-11.  Price  and  Plant  Mix  Under 
Diverse  Technology 

Note:  3-j  = capacity  cost  per  day 

32  = capacity  cost  per  day 
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peak-period  efficient  prices  and  higher  off-peak  efficient 
prices. 

(ii)  optimal  plant  mix  (x^  = off-peak  demand,  x 2 = peak  demand, 

= Plant  1 output,  and  q2  = Plant  2 output). 

The  cost  conditions  for  m technologies  are  as  follows: 

$1  ■>  32  ’ ■ ■ > ^ 5 

0 < b,  < . . . < b . 

1 m 

Capacity  is  installed  and  operated  in  order  of  increasing  operating 
cost.  In  particular,  Plant  I is  used  in  every  period  to  meet  all  de- 

mand. If  m = 2,  Plant  2 is  filling  the  remaining  demand.  Thus, 

= x^,  and  q2  = x2  - x,. 

To  illustrate,  assume  a plant  mix  consisting  of  two  plants  with  per 
unit  operating  costs  (b^  and  b^)  and  daily  per  unit  capacity  costs  (3j 
and  32) , 
where 

bj.  = 1.4,  b^  = 2,  3^  = 3,  and  32  = 2. 

The  demand  functions  are 

Pl  = ioo‘-xl/4> 

P2  = 200  -10X2 . 

The  welfare  optimal  prices  are 

px  = (2b1  + 3j_)  - (b2  + 32)  = 1-8,  and 
P2  _ (t>2  82 ) ~ ^ • 

Then  we  can  derive  x^,  x2  from  the  demand  function:  x^  = 16.07  and 

x2  = 19.60. 

It  is  straightforward  to  show  that  the  optimal  amounts  of  each 
plant  (q^  and  q2)  to  construct  are 
q^  = x^  and  q2  = x2  - x^, 
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with  Plant  1 being  operated  to  full  capacity  in  both  periods  and  Plant  2 
being  used  only  in  the  peak  period. 

In  Figure  3-11,  there  are  several  shaded  areas,  each  with  specific 
economic  interpretations: 

(i)  ATSV  represents  net  improvement  in  peak  period  allocative 
efficiency  (capacity  savings  which  outweigh  customer 
valuation  of  the  lost  output), 

(ii)  AABC  represents  net  consumer  surplus  gain  (after  taking 
the  transfer  into  account), 

(iii)  Peak  price  is  4<£/KWH;  however,  the  cost  of  peak  elec- 
tricity production  by  Plant  1 is  4.4<f/KWH.  The  total  loss 
is  represented  by  p?HJK.  In  contrast,  off-peak  price  is 
1.8d/KWH.  Up  to  an  output  of  X,,  the  off-peak  cost  is 
1.4<£/KWH;  FEBG  represents  total  revenues  greater  than 
costs  (quasi-rents  earned  from  this  consumption). 

This  model  of  diverse  technology  is  used  in  the  next  chapter  to 
generate  cost  data  relating  production  and  reliability.  The  concepts 
developed  here  establish  the  basis  for  formulating  a process  model  which 
permits  us  to  generate  pseudo  data  to  be  used  in  an  econometric  analysis 
of  costs. 


CHAPTER  4 

MODELING  THROUGH  A SIMULATION 'APPROACH 


This  study  differs  from  previous  studies  in  three  respects.  The 
first  is  that  the  reliability  of  power  supply  becomes  an  essential 
variable  in  the  analysis  of  electricity  costs.  In  most  previous  eco- 
nomic literature,  reserve  margin  and  the  probability  that  demand  does 
not  exceed  capacity  are  the  only  reliability  indices  used.  However,  two 
other  reliability  indices  widely  employed  by  engineers.  Loss  of  Load 
Probability  (LOLP)  and  Expected  Demand  not  Served  e(DNS),  are  considered 
here. 

The  second  contribution  we  will  attempt  to  make  involves  treating 
both  demand  and  supply  as  stochastic  variables.  Most  previous  studies 
consider  demand  as  the  only  stochastic  element,  which  severely  limits 
analysis,  especially  in  the  context  of  unit  outages.  The  third  aspect 
of  the  study  warranting  attention  is  the  use  of  pseudo-data  (the  results 
of  a set  of  simulations).  The  numbers  allow  us  to  apply  econometric 
techniques  to  determine  the  cost  of  reliability. 

4.1  The  Simulation  Model  (Process  Model) 

The  result  of  this  simulation  test  will  facilitate  several  economic 
analyses.  First,  the  data  permit  an  econometric  estimation  of  a cost  of 
reliability  equation: 


C - f (R,  TKWH,  K,  u), 
R = g(pKWH,  K,  v) , 


(4-1) 

(4-2) 
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where 
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C = total  cost, 

R = reliability, 

TKWH  = total  electricity  consumption, 
pKWH  = peak  load  electricity  consumption, 

K = capacity  generation  mix,  and 
u,  v = disturbance  terms. 

In  addition,  the  data  allow  us  to  determine  a set  of  prices  which  will 
maximize  the  net  benefit  subject  to  reaching  a given  level  of  relia- 
bility (measured  in  terms  of  a particular  index). 

4.1.1  Variables  and  Parameters  in  the  Model 

Before  getting  into  the  detailed  mechanism  of  the  simulation  model, 
a flow  chart.  Figure  4-1,  is  presented  to  summarize  the  interrelation 
among  relevant  variables  and  parameters.  This  flowchart  will  enable 
readers  to  follow  the  general  framework  of  this  simulation. 

Note  that  the  analysis  can  be  divided  into  two  areas.  The  first 
area  involves  a simulation  model  (process  model)  which  is  designed  to 
generate  pseudo  data  based  on  underlying  economic  and  engineering  rela- 
tionships. These  pseudo  data  are  generated  by  a well-specified  model. 
Then  we  use  alternative  regression  models  to  estimate  the  cost  of  relia- 
bility. The  econometric  tests  will  be  presented  in  the  next  chapter. 
For  now,  we  identify  the  underlying  assumptions.  Assume  that  demand  has 
the  following  functional  form: 

Di(pi , u.)  = f(pi)  + ui  (4-3) 

where 

p = price,  and 


u = disturbance  term. 
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Figure  4-1.  Flowchart  of  Simulation  Model 
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Demand  is  subject  to  random  fluctuations  so  it  is  stochastic.  The 
utility  company  (producer)  cannot  choose  a particular  price  and  be 

certain  as  to  what  the  load  response  (demand)  will  be.  We  assume  refer- 

ence demands  for  the  peak  and  off-peak  periods.  Instead  of  the  riskless 
relationship  D = D(p),  the  electricity  producer  is  confronted  with  D = 
D(p,  u).  Of  course,  demand  for  electricity  can  be  a function  of  elec- 
tricity price  (base  rate  and  fuel  adjustment),  per  capita  income,  popu- 
lation, and  other  factors.  For  simplicity,  we  assume  price  is  the  only 
independent  variable.  Also,  both  peak  and  off-peak  demands  are  in 
linear  form: 

D*(p*,  u*)  = a*  + d*p*  + u*  (4-3a) 

D 1 ( p 1 , u1)  = a'  + d'p'  + u1  (4-3b) 

where 

D*  = peak  demand  = q*, 

D'  = off-peak  demand  = q", 
p*  = peak  price, 
p1  = off-peak  price,  and 
u*,  u‘  = disturbance  terms. 

We  let  total  cost  be  a separable  function  of  capacity  cost,  and  of 
fuel  and  maintenance  cost.  These  two  costs  per  unit  of  output  are 
defined  as  3 and  b,  respectively. 

This  type  of  cost  function  has  been  used  extensively  in  the  litera- 
ture on  peak-load  pricing*.  We  further  assume  that  the  generation 


*The  average . capaci ty  cost  per  unit  of  output  assumes  that  the 
plant  is  fully  utilized.  Brown,  G.,  and  Johnson,  B.  M. , American  Eco- 
nomic  Review,  Vol . 59  (March  1969),  pp.  119-128. 
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capacity  is  divisible  for  running  purposes.  Also,  this  is  a closed 
utility  model,  so  we  do  not  consider  power  pooling.  For  simplicity, 
this  simulation  only  considers  fuel  cost  and  capacity.  We  also  assume 
that  the  base  unit  is  used  to  meet  demand  during  the  peak  period,  and 
the  peak  unit  is  used  to  fill  the  remaining  peak  demand.  Accordingly, 
the  total  cost  equations  for  firm  peak  and  off-peak  demands  under  two 
technologies  can  be  formulated  as  follows: 

C*  = (b‘  +3')q?  + (b*  + 3*)q*  = peak  total  cost,  (4-4a) 

U P 

C'  = b 1 q ' = off-peak  total  cost,  (4-4b) 

where 

C*  = peak  total  cost, 

C1  = off-peak  total  cost, 

3*  = peak-unit  capacity  cost  ($/KWH), 

3'  = base-unit  capacity  cost  ($/KWH), 
b*  = peak  fuel  cost  ( $/ KWH ) , 
b'  = base-load  fuel  cost  ($/KWH), 

-k 

qb  = peak  demand  supplied  by  base  unit  under  peak  load  pricing, 
q*  = peak  demand  supplied  by  peak  unit  under  peak  load  pricing, 
q*  = peak  demand  =q£  + q*,  and 
q1  = off-peak  demand. 

The  following  social  welfare  function  can  now  be  employed.  The 
expected  value  of  welfare  is 

E(W)  = E(TR  + CS  - TC),  (4-5) 

where 

TR  = total  revenue, 

CS  = consumer  surplus,  and 


TC  = total  cost. 
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By  Equation  4-3,  our  demand  function  in  general  form  is 

q = D(P,  y).  In  general,  the  demand  curve  can  also  be  written  as 

P = ^(q,  e).  Thus,  E(W)  = E(TR  + CS  - TC)  can  be  written  as 

E(W)  = (p*q*  - (/jj*  (q*,  e*)dq*  - p*q*)  - C*) 

+ ( p 1 q ' - (/  jj*(q',  e')dq'  - p'q')  - C1)  (4-5a) 

where 

p*  = peak  price, 
p'  = off-peak  price,  and 
y,  e = disturbance  terms. 

A practical  question  can  now  be  posed:  would  it  be  socially  effi- 

cient to  switch  from  uniform  pricing  to  time-of-day  pricing?  Our  gener- 
al approach  to  this  problem  may  be  treated  as  the  problem  of  maximizing 


a net  welfare  function  of  the  following  form: 

Wx  = TRj  + CS1  - TCj,  (4-6) 

W2  = TR2  + CS£  - TC2,  (4-7) 

W2  - W1  = (TR2  - TRX)  - (TC2  - TC1  ) + (CS2  - CS^,  (4-8) 

AW  = Ait  + ACS,  (4-9) 

TC1  = C**  + C‘ 1 , (4-6a) 

TC2  = C*  + C' , (4-7  a ) 


where  subscript  1 represents  the  uniform  pricing  case,  and  subscript  2 
represents  the  time-of-day  pricing  case, 
c**  = (b1  + 3')q£*  + (b*  + 8*)q**> 
c"  = b'q"  , 

q**  = peak  demand  supplied  by  base  unit  under  uniform  pricing, 
q**  = peak  demand  supplied  by  peak  unit  under  uniform  pricing,  and 


q11  = off-peak  demand  under  uniform  pricing. 
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The  term  AW  represents  net  welfare  returns  (net  benefit)  under  perfect 
rationing  (costless  rationing  according  to  willingness  to  pay).  In 
addition,  Air  represents  the  profit  change  between  uniform  price  and 
time-of-day  prices  (i.e.,  peak  price  and  off-peak  price).  The  term  ACS 
consists  of  the  net  consumer  surplus  changes.  One  change  is  from  the 
difference  between  peak  price  and  base  price  (uniform  price)  while  the 
other  is  from  the  difference  between  off-peak  price  and  base  price. 

We  can  also  incorporate  a reliability  constraint  into  the  above 
analytical  framework.  In  addition,  the  derivation  of  an  optimal  set  of 
time-of-day  prices  can  be  obtained  through  the  maximization  process  of 
the  above  welfare  function,  subject  to  reaching  a given  level  of  relia- 
bi  1 ity. 

The  general  form  of  the  reliability  equation  can  be  written  as 
follows: 

R = g(D*,  K)  + v,  (4-10) 

where 

R = reliability  index  (LOLP  or  e(DNS)  or  RM), 

K = capacity  generation  mix, 

D*  = peak  load,  and 

v = disturbance  term. 

Capacity  generation  mix  refers  to  different  combinations  of  genera- 
tors. Each  of  the  generators  has  its  own  capability,  fuel  type,  and 
outage  rate.  Thus,  different  capacity  mixes  imply  different  generation 
reliability.  To  calculate  LOLP  or  e(DNS),  two  items  of  information  are 
required: 

(a)  units  in  the  generating  system  (supply  side):  capability,  and 

forced  outage  rate,  and 
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(b)  the  pattern  of  demand  (demand  side).  (For  simplicity,  we  will 
summarize  the  pattern  in  the  simulation  with  one  parameter, 
the  daily  peak  load.) 

In  summary,  our  simulation  model  consists  of  Equations  4-3a,  4-3b, 
4-4a,  4-4b,  4-5a,  4-9,  and  4-10.  The  data  simulation  framework  is  built 
on  this  system  of  seven  equations.  We  discuss  the  details  of  the  simu- 
lation procedure  and  data  base  in  the  following  section. 

4.1.2  Data  and  Simulation  Procedures 

The  rating  periods  used  to  define  peak  and  off-peak  in  this  simula- 
tion test  are  illustrated  below: 

12  midnight  ) off-peak 

8 a,m-  j peak  (10  hrs. ) 

6p.m.  | 

12  midnight  ) off-peak 

Thus,  the  reference  load  pattern  is  very  simple,  as  shown  in  Figure  4-2. 
The  dotted  lines  represent  the  net  result  of  two  factors:  price  struc- 

ture and  random  demand  shocks. 

To  obtain  peak  and  off-peak  load  patterns,  we  use  a random  number 
generator  to  provide  normally  distributed  load  data  for  a reference  load 
shape:  peak  load  (with  a mean  of  instantaneous  peak  load  Xp  = 46.5MW 

and  a standard  deviation  Sp  = 5.5  MW)  and  off-peak  load  (mean  xQp  = 23.5 
MW  and  a standard  deviation  sQp  = 5.5  MW).  Thus,  each  day  (for  a 
thirty-day  period)  has  a different  load  shape.  To  achieve  annualized 
numbers,  the  monthly  KWH  observation  was  multiplied  by  twelve.  The  peak 
instantaneous  KW  demand  was  also  used  in  the  regression  analysis.  Thus, 
long  run  capacity  adjustments  and  seasonal  variations  in  generating  unit 
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MW 


8am  6pm  12pm 


Hour 


Figure  4-2.  Reference  Load  Patterns  and  Shifts  Due 
to  Prices  and  Random  Shocks 
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maintenance  or  load  patterns  are  not  being  considered  here.  The  same 
reference  load  pattern  is  used  for  each  generation  mix.  The  random 
loads  were  created  by  using  the  reference  load  adjusted  by  a random 
factor  multiplier  <f>  = (e)0.00005(  (-1  )mn)  (mn) , 
where 

tf>  = the  random  multiplier  determining  the  size  of  the  shift, 
e = the  elasticity  adjustment  factor, 

■f*  h 

m = the  nr  price  combination  (m  = 1 . . . 60),  and 
n = the  n^  day  (n  = 1 . . . 30). 

Thus,  the  reference  load  is  shifted  up  or  down,  depending  on  the  day  and 
price  combination.  The  size  of  the  shift  (relative  to  the  reference 
load)  depends  on  the  price  combination,  the  associated  elasticities,  and 
a purely  random  factor. 

We  have  developed  two  major  groups  of  scenarios  which  vary  in  price 
and  demand  elasticity  (own  and  cross).  All  relevant  parameters  of  two 
major  groups  mentioned  above  are  summarized  in  Table  4-1.  (The  deriva- 
tion of  these  parameters  is  fully  discussed  throughout  the  rest  of  this 
chapter.)  Each  major  group  is  composed  of  four  subgroups  based  on 
different  capacity  mixes,  which  are  specified  in  Figure  4-3. 

The  decision  regarding  the  selection  of  unit  size  and  numbers  of 
units  is  based  on  economic  and  engineering  considerations.  The  former 
refers  to  load  patterns  and  load  growth,  capacity,  and  fuel  costs.  The 
latter  refers  to  engineering  efficiency.  (The  full  treatment  of  this 
issue  is  beyond  the  scope  of  this  study.)  For  this  simulation,  the 
selection  of  reference  capacity  mix  (x-mix)  is  based  on  the  following 
load  considerations.  As  we  generated  peak  and  off-peak  demand  in  the 


Table  4-1.  Creation  of  Pseudo  Data 
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(I)  Reference  Mix 


4 (cost/kwh) 


Bl 

b2 

A 

22 

30 

18 

0 22  52  70 


Capacity  (MW) 


(2)  Equal  Size  Base  Units 


4 (cost/kwh) 


Capacity  (MW) 


(3)  Relatively  More  Peaking  Capacity 


4(cost/kwh) 


1.169 

0.928 
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B2 
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Capacity  (MW) 


(4)  High  Total  Capacity 
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Figure  4-3.  Four  Generating  Mixes 

(Minimal  Capacity  Variance) 
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earlier  paragraphs,  we  have 

(i)  Xp  = 46.5  MW  sp  = 5.5  MW; 

therefore 

p ( x < x + s = 46.5  MW  + 5.5  MW  = 52MW)  = 84%. 

P P P 

We  make  total  base  capacity  equal  to  the  mean  plus  one  standard  devia- 
tion (i.e.,  + B2  = 52  MW).  We  assume  that  Unit  = 22  MW,  Unit  B2  = 

30  MW. 

(ii)  xQp  = 23.5  MW  sQp  =5.5  MW, 

therefore 

p(xop  < xQp  + 3sop  =”23.5  MW  + 3(5.5  MW)  = 40  MW)  = 99%  . 

We  have  to  make  Unit  A = 18  MW  in  order  to  satisfy  off-peak  demand  if 
either  B^  or  B2  is  out  because  of  forced  outage  or  scheduled  mainte- 
nance. The  following  four  types  of  capacity  variations  (see  Figure  4-3) 
were  created  to  explore  the  cost  implications  of  alternative  capacity 
mixes. 

1.  Reference  k-mix:  total  capacity  = 70  MW, 

Base  Units:  Unit  B.  (22  MW),  Unit  B?  (30  MW),  and 

Peak  Unit:  Unit  A t 18  MW); 

2.  Equal  size  Base  Units:  total  capacity  = 70  MW, 

Base  Units:  Unit  B,  (26  MW),  Unit  B?  (26  MW),  and 

Peak  Unit:  Unit  Ai(18  MW); 

3.  Mix  of  Peak  Unit  and  Off-Peak  Unit:  total  capacity  = 70MW, 

Base  Units:  Unit  B.  (18  MW),  Unit  (30  MW),  and 

Peak  Unit:  Unit  Ai(22  MW)  (i.e.,  relatively  more  peaking 

capacity) ; 

4.  Change  in  Total  Capacity:  total  capacity  = 77  MW  (i.e., 

10  percent  capacity  increase  for  all  units), 

Base  Units:  Unit  B-,  (24.2  MW),  Unit  B?  (33  MW),  and 

Peak  Unit:  Unit  A1(19.8  MW). 


The  base  price  (uniform  price)  is  assumed  to  be  6 <t  per  KWH,  which 
happens  to  be  the  average  electricity  price  in  the  State  of  Florida  in 
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1981.  In  the  simulation,  60  different  price  combinations  were  used  for 
each  generation  mix.  Each  price  combination  represents  one  set  of  peak 
and  off-peak  price.  The  price  differences  between  peak  and  off-peak 

price  are  specified  as  follows: 

1.  For  the  first  20  price  combinations  the  peak  and  off-peak  price 
differences  are  1.5<£, 

2.  For  the  second  20  price  combinations  the  peak  and  off-peak  price 
differences  are  2<t,  and 

3.  For  the  last  20  price  combinations  the  peak  and  off-peak  price 
differences  are  4<£. 

The  set  of  time-of-day  price  elasticities  were  selected  from  a load 
study  by  Cargill  and  Meyer  (1971).  The  data  base  of  Cargill's  study 
consists  of  monthly  observations  on  the  twenty-four  hour  load  curve  in 
one-hour  intervals  for  two  regions  during  the  four-year  period  from 
January  1965  through  December  1968.  The  first  region  is  a large  midwest 
industrial  area  and  the  second  is  a west  coast  agricultural  region. 
Their  estimated  price  elasticities  for  each  hour  are  presented  in  Table 
4-2. 

Since  our  peak  hours  are  defined  from  8 a.m.  to  6 p.m.,  the  peak 
and  off-peak  elasticities  are  the  average  hourly  elasticities  in  two 
periods  (peak  and  off-peak).  These  two  elasticities  are  taken  to  be  as 
fol 1 ows : 

Region  1 Region  2 

Peak  0.492  0.241 


Off-peak 


0.487 


0.315 
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Table  4-2.  Time  of  Day  Pri'ce  Elasticities 


Hour 

Region  1 

Region  2 

1 a .m. 

-0.50 

-0.49 

2 

-0.50 

-0.52 

3 

-0.51 

-0.42 

4 

-0.51 

-0.48 

5 

-0.50 

-0.48 

6 

-0.49 

-0.29 

7 

-0.58 

-0.27 

8 

-0.56 

-0.23 

9 

-0.56 

-0.27 

10 

-0.57 

-0.35 

11 

-0.53 

-0.34 

12  noon 

-0.48 

-0.26 

1 p.m. 

-0.43 

-0.19 

2 

-0.47 

-0.20 

3 

-0.44 

-0.20 

4 

-0.43 

-0.18 

5 

-0.46 

-0.09 

6 

-0.44 

-0.18 

7 

-0.45 

-0.06 

8 

-0.57 

-0.15 

9 

-0.55 

-0.13 

10 

-0.45 

-0.23 

11 

-0.39 

-0.30 

1 2 midnight 

-0.38 

-0.35 

Data  Source: 

Cargill  and  Meyer,  (1971). 
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Cross-elasticities  used  in  these  scenarios  are  presented  as 
follows: 


Peak 

Off-peak 


Region  1 Region  2 

0.01  0.02 

0.06  0.05 


1 

\ 

4.1.3a  Welfare  Calculations 

Now  we  turn  to  the  discussion  of  estimation  of  consumer  surplus. 
In  Equation  4-9,  net  change  in  consumer  surplus  caused  by  price  change 
can  be  represented  by  the  area  P P 1 1 GE  in  Figure  4-4,  which  can  be 
approximated  by  the  following  equation: 

ACS  = ( P - P ' ' )Q  + %((P  - P")(Q' 1 - Q)) 

- (P  - P")(Q  + %(Q"  - Q)) 

- h(P  - P")(Q  + Q"),  (4-11) 


where 


ACS  = net  change  in  consumer  surplus, 

P = price  before  price  change, 

Q = quantity  before  price  change, 

P11  = price  after  price  change,  and 
Q'1  = quantity  after  price  change. 

The  above  estimation  is  based  upon  the  assumption  that  the  demand 
function  is  linear  and  independent. 


4.1.3b  Capacity  Cost  Calculations 

The  selection  of  generator  unit  size  (i.e.,  capacity  level)  and 
number  of  generation  units  are  determined  by  historical  load  patterns. 
The  simulation  results  reported  here  do  not  reflect  the  cost  of  optimal 


Demand  Function 


F 


0 Q Q" 


KWH 


Figure  4-4.  Net  Change  in  Consumer  Surplus 
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generating  mixes  for  the  ex  post  load  pattern.  This  characteristic  is 
not  a severe  limitation,  since  the  generating  mixes  of  most  electric 
utilities  were  non-optima!  in  the  late  1970s,  given  the  input  prices 
after  the  OPEC-induced  petroleum  price  rises.  The  stochastic  demands 
are  still  met  with  minimum  operating  costs.  So  the  costs  and  outputs 
observed  are  not  unlike  the  data  actually  reported  by  electric  utilities 
during  similar  time  periods. 

Since  we  are  dealing  with  an  extremely  simple  production  technol- 
ogy, we  further  assume  that  the  generator  capacity  is  divisible.  This 
simplifying  assumption  does  not  violate  the  essential  feature  of  a 
diverse  technology  production  environment.  Also,  note  that  this  utility 
model  assumes  a closed  system,  so  it  does  not  consider  electricity 
imports  or  exports.  That  is,  power  pooling  issues  do  not  arise. 

In  this  simulation,  we  assume  there  are  three  generating  units. 
Unit  and  Unit  B2  are  coal-fired  base  units.  Unit  A is  the  oil-fired 
peaking  unit.  The  capacity  cost  can  be  defined  in  terms  of  the  concepts 
developed  in  the  previous  chapter. 

The  annual  capacity  cost  if  unit  i,  denoted  by  3*,  is  defined  as 
follows : 

= FicciKi,  (4-12) 

where 

F.j  = fixed  charge  rate, 

cc.  = generator  capacity  acquisition  cost  ($/MW), 

K.j  = capacity  in  MW,  and 

3*  = total  annualized  capacity  cost. 

We  use  the  above  definition  to  calculate  an  annual  capacity  cost. 
Further,  if  Equation  4-12  were  divided  by  the  expected  value  of  the 


J 
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annual  energy  produced  by  unit  i,  denoted  by  e(E^),  then  the  capacity 
cost  might  be  expressed  in  Mils/ KWH  or  $/MWH. 

We  have 


8*  F.cc.k.  F.cc.k. 

HTJ  = ^lETT  = 3,650(CF.k.)  > <since  we  have  10  hours  peak  in 
1 1 IT 

simulation,  so  e(E^)  - 3,650(CF^k_j ) ) , 


pi  = 


Ficci 

3,650(CFi) 


(4-13) 


where 

CF.j  = capacity  factor. 

The  capacity  factor  is  the  ratio  of  the  average  load  to  the  total 
plant  capacity.  Given  the  simultaneity  problem  linking  KWH  consumption 
to  this  number  we  will  assume  that  CF  = 0.8  for  all  base  units,  CF  = 0.2 
for  the  peaking  unit,  and  the  generator  acquisition  costs  (unit  capacity 
costs)  are  assumed  as  follows: 


UC  (30  MW)  = $190/KW, 

UC  (18  MW)  = $1 95/KW , and 
UC  (22  MW)  = $ 65/KW . 

Then,  we  use  Equation  4-13  to  calculate  capacity  cost  (8^)  as  follows:* 
(i ) Base  Uni ts: 


Ficci 


30MW 


3,650(cf  ry 


0. 1 79(1 90)1 03 
3,650(0.8) 


= 11.49  Mi  11  s/KWH  = 1.149<£/KWH, 


*We  make  the  following  assumptions  regarding  the  generator  units: 
1)  straight-line  depreciation,  2)  life  of  generator  is  assumed  to  be  10 
years,  and  3)  zero  net  salvage  value.  Also,  we  assume  the  depreciation 
rate  = 10%,  the  tax  rate  = 5.7%,  other  local  and  state  taxes  = 1.6%,  and 
the  insurance  rate  = 0.6%.  Thus,  the  sum  of  the  above  percentages  is 
17.9  percent,  and  the  fixed  charge  rate  (FCR)  is  0.179  for  all  units. 


'■nJ 
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where  F..  = 0.179; 

„ _ 0.179(195)103 

H8MW  3,650(0.8) 


1 .188(t/KWH. 


Thus,  a base  unit  capacity  cost  of  1.169<£/KWH  has  been  used  in  this 
simulation  test. 

( ii ) Peaki ng  Uni t: 

a - 0- 1 79(65)1 03  _ n QOQQj./1/im 
322MW  3,650(0.2)  0 . 9288<t/KWH . 


Fuel  prices  are  major  factors  in  determining  the  production  cost. 
According  to  the  1980  fuel  price  level,  it  costs  1.877  <t/KWH  for  a coal 
unit  (10,485  BTUs  x $1.8  per  million  BTU  = 1.887  tf/KWH),  and  3.902  <f/KWH 
for  an  oil  unit  (10,485  BTUs  x $3,722  per  million  BTU  = 3.902  <£/KWH). 
According  to  composite  averages  (all  fuel  consumed  for  the  State  of 
Florida  in  the  year  1980),  BTU  per  net  KWH  is  equal  to  10,485.* 

Three  reliability  indices  are  simultaneously  used  in  this  simula- 
tion test:  Expected  demand  not  served  e(DNS),  loss  of  load  probability 
(L0LP),  and  reserve  margin  (RM).  Since  load  may  be  lost  due  to  many 
causes,  and  only  a shortage  of  generating  capacity  is  considered  in  the 
L0LP  method.  This  index  is  measured  in  expected  days  of  capacity  short- 
age per  year,  and  is  often  referred  to  as  a system's  risk  index.  The 
risk  index  (L0LP),  is  affected  by  the  installed  capability  and  peak  load 


*The  British  Thermal  Unit  (BTU)  is  the  standard  unit  for  measuring 
the  quantity  of  heat  energy,  such  as  the  heat  content  of  fuel.  It  is 
defined  as  the  amount  of  heat  energy  necessary  to  raise  the  temperature 
of  one  pound  of  water  one  degree  Farenheit. 
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just  as  is  the  percent  reserve  index;  however,  the  LOLP  is  also  sensi- 
tive to  the  various  sizes  and  outage  rates  of  the  installed  generator, 
the  load  shape,  and  maintenance  scheduling.  For  this  simulation  test, 
the  following  table  provides  the  necessary  information  to  calculate  loss 
of  load  probability. 

Table  4-3  contains  the  forced  outage  parameters  used  in  this  simu- 
lation. The  largest  unit,  E^,  is  assumed  to  have  the  highest  forced 
outage  rate  where  the  definition  of  forced  outage  rate  is 

Forced  _ forced  unavailable  hours 1 

outage  rate  service  hours  + forced  unavailable  hours  ' ' 

In  a more  realistic  modeling  exercise,  unit  forced  outage  rates 
would  be  taken  from  the  referenced  EEI  publication  or  alternatively  from 
a utility  company's  historical  operating  records  on  the  unit  in  ques- 
tion. We  assume  independence  of  unit  outage,  so  a cumulative  outage 
probability  table  can  be  calculated.  (See  Table  4-4  for  the  appropriate 
figures. ) 

To  calculate  reliability  indices,  the  generating  system  is  repre- 
sented by  a capacity  outage  table,  which  requires  enumerating  all  pos- 
sible combinations  of  the  units  in  or  out  of  service  and  the  associated 
probability  of  occurrence  for  each  combination.  Table  4-4  contains  an 
outage  probability  table  which  shows  both  exact  outage  (the  particular 
outage  case)  and  cumulative  outage  probability. 

As  to  the  calculation  of  expected  demand  not  served,  e(DNS)  can  be 
defined  as  follows: 

e(DNS)  = E(D  - S) 

= E(x) 

= 5 x,  f(x,), 

i=l  1 1 


(4-15) 
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Table  4-3.  Forced  Outages  for 
Units  in  Generating  System 


Uni  t 


A 

B, 


Capability  (MW) 


Forced  Outage  Rate 


22 


0.01 


18 


0.02 


Inage  Rate 
0.99 
0.98 


30 


0.03 


0.97 


Table  4-4.  Exact  Outage  and  Cumulative  Outage  Probability  Table 
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where 

D = peak  load , and 
S = supply  of  electricity. 

Both  D and  S are  random  variables,  and  if  we  let  x = D - S,  x 
becomes  a discrete  random  variable.  As  an  illustration,  consider  the 
following  example. 

Let  the  peak  load  on  the  7th  day  be  21  MW.  What  is  the  expected 
demand  not  served?  An  outage  of  more  than  49  MW  (i.e.,  total  available 
capacity  (70  MW)  minus  peak  load  (21  MW)  equals  49  MW)  will  cause  a loss 
of  load.  In  other  words,  only  Cases  7 and  8 in  Table  4-4  will  cause  a 
loss  of  load.  Therefore,  if  Case  7 occurs,  then  DNS  = 21  - 18  = 3 MW, 
where 

DNS  = peak  load  - capacity  in  service. 

A1  so , 

f(x})  = Prob  (Case  7)  = 0.000294. 

If  Case  8 occurs , then 

DNS  = 21  - 0 = 21  MW; 

f(x£)  = Prob  (Case  8)  = 0.000006. 

By  using  Equation  4-15  we  obtain 
2 

E(X)  = E X.  f ( x • ) = 3 x 0.000294  + 21  x 0.000006  = 1.008  KW. 
i=l  1 

Accordingly,  the  expected  demand  not  served  for  the  7th  day  is  1.008  KW 
with  a probability  of  0.0003.  Since  both  e(DNS)  and  L0LP  are  created 
from  the  same  set  of  data,  we  are  now  able  to  understand  that  one  side 
of  the  two-way  mirror  is  L0LP;  the  other  side  is  e(DNS),  measured  in  KW 
per  time  period. 

The  summary  of  all  the  relevant  parameters  used  in  this  simulation 
test  is  found  in  Table  4-5.  Relevant  parameters  for  two  groups  of 


Table  4-5.  Parameter  List  of  Simulation  Test 
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scenarios  were  listed  separately.  These  parameters  are  uniform  price, 
price  elasticity,  price  cross  elasticity,  capacity  charge,  and  fuel 
cost. 

Up  to  this  point,  we  have  specified  and  described  all  the  required 
simulation  procedures  and  data  base.  All  the  simulation  results  and 
further  econometric  analysis  are  presented  in  the  next  chapter. 


CHAPTER  V 

EMPIRICAL  ESTIMATES  OF  THE  ELECTRICITY 
COST  OF  RELIABILITY  AND  PRODUCTION 

We  now  review  the  potential  significance  of  the  simulation  results 
by  performing  some  fundamental  econometric  analyses  and  making  numerical 
comparison.  We  are  mainly  concerned  with  the  estimation  of  cost  of 
reliability  by  using  pseudo  data.  We  are  also  interested  in  some  impli- 
cations of  the  simulation  for  non-cost  dimensions  of  performance.  For 
example,  we  can  obtain  a set  of  welfare-maximized  TOU  prices  subject  to 
reaching  a given  level  of  reliability  constraint. 

5,1  Structure  of  the  Model 

The  first  major  task  in  this  section  is  to  conduct  an  econometric 
analysis  based  on  results  from  the  simulation  test.  That  is,  we  need  to 
estimate  the  cost  of  reliability  equations  to  assist  us  in  power  plan- 
ning and  investment  decisions.  By  examining  Table  4-4,  we  see  that 
there  are  eight  possible  outage  combinations.  Each  of  the  situations 
involves  different  sets  of  generating  unit  status  (i.e.,  up  or  down). 
Therefore,  each  of  the  cases  has  a different  combination  of  production 
costs  and  reliability  for  a given  load  pattern. 

5.1.1  Determinants  of  Reliability 

Before  going  into  the  estimation  model,  we  will  give  a brief  dis- 
cussion of  determinants  of  reliability.  Most  power  outages  occur  during 
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the  peak  hour,  so  the  instantaneous  peak  KW  demand  will  affect  the 
reliability  of  an  entire  power  system.  The  reliability  is  negatively 
related  to  peak  electricity  demand  which  is  also  negatively  related  to 
peak  price.  Thus,  reliability  is  positively  related  to  peak  price. 

Load  factor  is  defined  as  the  ratio  of  the  average  load  in  KW 
supplied  during  a designated  period  to  the  peak  or  maximum  load  in  KW 
occurring  in  that  period.  Accordingly,  the  higher  the  load  factor,  the 
lower  the  system  reliability  will  tend  to  be,  for  the  same  generation 
mix. 

i 

5.1.2  The  Model 

The  functional  relationship  of  the  model  is  formally  described  as 
follows: 

C - f (X,  R,  e),  and  (5-la) 

R = g(p,  L,  p),  , (5-2a) 

where  C denotes  the  dependent  variable,  total  cost,  and 
X = total  electricity  consumption, 

R = reliability, 
p = price  of  electricity, 

L = peak  load,  and 

e,  p = disturbance  terms. 

Regressor  R in  the  above  system  of  equations  is  endogenous,  sto- 
chastic, and  correlated  with  the  regression  disturbance.  In  this  case, 


'Sj 
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the  ordinary  least  squares  estimators  of  the  regression  coefficients  are 
inconsistent.*  There  are  several  methods  leading  to  consistent  estima- 
tion. For  simultaneous  equation  systems  (including  the  recursive  sys- 
tem), the  best  known  common  method  is  that  of  two-stage  least  squares. 
In  considering  Equation  5-2a,  L and  P are  closely  correlated**  because 
L-j  = l(P.j);  we  need  to  drop  P so  that  Equation  5-2a  is  modified  for 
estimation  purposes  as  follows: 

Ri  = Y0  + n ■ Lijk  + »i-  ( 5— 2b) 
We  can  also  rewrite  Equation  5-la  in  this  form: 


Ci  = e0  + 61  • Xijk  + b2  - Rijk  + ei 


(5-lb) 


where 


k-mix 


= 1 ...  4, 

j = elasticity  combination 
= 1,  2,  and 
k = price  combination 
= 1,  60. 


On  a priori  grounds,  we  may  hypothesize  the  sign  for  some  of  the 
variables.  In  the  cost  equation,  we  expect  positive  marginal  production 


*See  Kmenta,  "Elements  of  Econometrics,"  Chapter  13.  In  Section 
13-1,  he  emphasized  that  the  ordinary  least  squares  method  of  estimation 
applied  to  the  structural  equations  of  a simultaneous  equation  system, 
and,  in  general,  leads  to  inconsistent  estimates.  A high  degree  of 
mul ticol 1 i neari ty  exists  whenever  one  explanatory  variable  is  highly 
correlated  with  another  explanatory  variable  or  with  a linear  combina- 
tion of  other  explanatory  variables.  Multicol 1 inearity  implies  that  the 
higher  the  degree  of  multicoil inearity,  the  larger  the  variances  of 
estimators  and  the  co-variances  between  estimators. 

**The  mul ti col  1 i neari ty  problem  is  not  precisely  eliminated  by 
suppressing  p,  rather,  y-1  will  be  a composite  coefficient  showing  the 
effects  of  both  p and  L. 
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cost,  so  the  sign  of  the  X coefficient  (total  electricity  consumption), 
31S  is  positive,  and  the  sign  of  the  R coefficient  (reliability),  32>  is 
expected  to  be  positive,  since  other  things  being  equal,  higher  reli- 
ability will  cause  higher  cost.  Therefore,  the  sign  of  R is  as  follows: 
if  R is  LOLP,  B2  < 0,  j 

if  R is  e(DNS),  B2  < 0,  and  ! 

if  R is  RM,  82  > 

In  the  reliability  equation,  the  sign  of  peak  load  is  hypothesized  as 
follows: 

if  R is  LOLP,  > 6, 

if  R is  e(DNS),  Yj  > 0,  and 

if  R is  RM,  < 0. 

Equations  5-lb  and  5-2b  belong  to  a system  of  linear  equations.  We 
can  test  whether  a nonlinear  model  will  yield  a better  fit.  Such  a 
system  of  non-linear  models  is  hypothesized  as  follows: 


We  can  also  transform  the  above  equations  into  a linear  model  in 
logarithmic  form  by  taking  logarithms  (to  base  e)  of  both  sides  of  Equa- 
tions 5-3  and  5-4: 


(5-3) 

(5-4) 


(5-5) 

(5-6) 


which  can  be  written  as 


(5-7) 

(5-8) 


•''NJ 
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where 


Ct  = In  C., 

x?jk  - '»  xijk- 

R?jk  ■ 1n  Rijk' 

Lijk  ■ Lijk- 
Bq  = In  8q,  and 


Y? 


In  yn. 


'0  'O' 

Here,  the  starred  symbol  represents  the  logarithms  of  the  unstarred 
counterparts. 


5.2  Econometric  Results 

We  apply  the  two-stage  least  square  (2SLS)  estimation  method  to 
both  the  linear  model  (i.e.,  Equations  5-lb  and  5-2b)  and  the  loga- 
rithmic model  (i.e.,  Equations  5-7  and  5-8).  Since  three  reliability 
indices  have  been  used  in  this  study,  we  will  estimate  the  equation  for 
each  index,  yielding  six  sets  of  estimations  for  each  capacity  mix 
(since  there  are  two  models,  linear  and  logarithmic,  and  three  relia- 
bility indices,  LOLP,  e(DNS),  and  RM). 


5.2.1  Hypothesis  Testing 

For  the  evaluation  of  econometric  results,  we  need  to  discuss 
methods  of  testing  the  hypotheses.  The  null  hypothesis  is  HQ:  8=0. 
We  wish  to  test  Hg  at  a one  percent  level  of  significance.  Since  we 
have  120  observations  for  each  K-mix,  the  appropriate  value  of  the  t- 
statistic  for  a two-tail  test  is  the  value  corresponding  to  t,,-,  n nnr, 
for  the  cost  equation,  and  t^g  g gg^  for  the  reliability  equation. 
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The  acceptance  region  for  our  test  is 


A 

-2.576  < < 2.576  . (5-9) 


Also,  the  tabulated  values  of  F(<_1  n_k  (i.e.,  F2  ^ and  F1  ^1Q)  at 
the  1 percent  level  of  significance  are  4.81  and  6.87  respectively. 

Based  on  these  results,  we  make  the  following  conclusions  regarding 

2 

R s,  t-values,  and  F-values: 

(1)  These  t-values  of  all  estimated  coefficients  of  both  linear 
and  logarithmic  models  are  out  of  the  acceptance  region  as 
indicated  by  Equation  5-9.  Thus,  all  estimated  coefficients 
are  statistically  significant  and  the  null  hypotheses  would  be 
rejected . 

(2)  The  F values  are  also  out  of  acceptance  region.  That  is,  the 
evidence  strongly  suggests  that  all  explanatory  variables 
significantly  contribute  to  the  explanation  of  variation  in 
cost. 

2 

(3)  All  R s are^high  enough  to  claim  a good  fit;  however,  there  is 
not  much  R^  variation  among  all  cost  equations,  that  is,  they 
are  all  over  99  percent,  suggesting  multicollinearity. 

(4)  The  logarithmic  model  has  a built-in  elasticity  measure.  For 
example,  in  Equation  5-7,  8,  measures  production  elasticity  of 
cost  and  measures  reliability  elasticity  of  cost.  We  will 
further  elaborate  such  measures  in  the  following  sections. 


5.2.2  Interpretation  of  Coefficients  on  Reliability  and  Cost 
We  now  turn  to  the  implications  of  estimated  models.  The  relia- 
bility elasticity  of  cost  is  defined  as 

r = 9(ln  C) 

^RC  9(ln  R)  ' 


Similarly,  the  electricity  production  cost  elasticity  is  defined  as 

r _ 9(ln  C) 

LXC  9 ( 1 n X)  * 


-"'J 
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We  can  derive  the  above  elasticities  from  a log  model.  Table  5-1 
presents  all  three  reliability  elasticities  of  cost.  The  meanings  of 
the  numbers  in  Table  5-1  are  very  straightforward.  For  instance, 
-0.0106  means  that  a 1 percent  decrease  in  LOLP  (i.e.,  reliability 
increases)  will  result  in  a 0.0106  percent  cost  increase.  (It  should  be 
noted  that  the  elasticities  in  Table  5-1  are  inelastic.) 

In  the  long  run,  an  improvement  of  reliability  can  be  obtained  by 

gC 

investing  in  additional  generating  units.  In  theory,  > 0 , 

(Both  cost  and  reliability  move  in  the  same  direction). 

In  the  short  run,  any  generator  outage  will  reduce  service  relia- 
bility and  total  electricity  production.  Since  the  capacity  cost  of  an 
out-of-service  unit  is  included  in  the  calculation  of  total  cost,  the 
total  cost  is  reduced  due  to  less  fuel  usage.  Thus,  > 0 , should 

still  hold  true  in  the  short  run.  In  Table  5-1,  the  signs  of  the  coef- 
ficients under  the  RM  column  are  all  negative.  However,  in  theory,  they 
should  be  positive.  This  discrepancy  suggests  that  if  the  reliability 
index  is  reserve  margin,  then  the  log  model  should  not  be  utilized.  We 
can  also  derive  the  following  findings  by  comparing  the  coefficients  in 
Table  5-1: 

(i)  For  all  four  types  of  k-mix,  costs  are  less  sensitive  to  the 
LOLP  reliability  index  than  to  e(DNS),  since  the  absolute 
value  at  each  LOLP  coefficient  is  smaller. 

(ii)  To  compare  k-mix  Type  1 (reference  k-mix)  and  k-mix  Type  2 
(two  equal  size  base  units),  the  latter  has  smaller  elas- 
ticity. That  is,  the  greater  the  generator  size  variation, 
the  greater  the  cost  sensitivity. 

(iii)  To  compare  k-mix  Type  1 (reference  mix)  and  k-mix  Type  3 

(larger  peaking  unit),  Type  3 k-mix  has  larger  elasticities, 
which  suggests  that  having  a larger  peaking  unit  in  the  capa- 
city mix  tends  to  cost  more  due  to  a given  reliability  change. 

(iv)  To  compare  k-mix  Type  1 (reference  mix)  and  k-mix  Type  4 

(10  percent  capacity  increase  for  all  units),  we  see  that  the 
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reliability  elasticity  of  cost  is  getting  larger  as  the  total 
capacity  of  generators  is  increasing. 

Table  5-2  summarizes  the  marginal  cost  of  production  (holding 
reliability  constant),  which  are  derived  from  our  linear  model.  We  can 
see  that  the  standard  errors  in  the  e(DNS)  column  are  the  smallest  among 
all  three  reliability  indices,  which  suggests  that  perhaps  the  e(DNS) 
model  gives  the  best  estimate  of  the  marginal  cost  of  electricity  pro- 
duction. By  comparing  k-mix  Type  1 with  k-mix  Type  4,  we  also  see  that 
marginal  costs  are  increasing  as  the  sizes  of  generators  are  increasing 
for  both  LOLP  and  e(DNS)  models. 

Table  5-3  summarizes  the  sensitivity  coefficients  of  reliability  to 
peak  demand,  that  is,  ||^— • These  coefficients  provide  insight  into 
reliability  change  with  respect  to  peak  demand.  Looking  at  the  ratio, 
we  can  see  that  the  e(DNS)  model  has  a substantially  lower  standard 
error  for  all  types  of  k-mixes  indicating  that  the  e(DNS)  model  may  pro- 
vide the  best  estimate  of  reliability  change  with  respect  to  peak  demand 
change.  Also,  by  comparing  k-mix  Type  1 with  k-mix  Type  4,  we  once 
again  see  that  k-mix  Type  4 has  larger  coefficients  for  both  LOLP  and 
e(DNS)  models.  This  means  that  reliability  of  k-mix  with  larger  units 
is  more  sensitive  to  peak  demand  change. 

Another  interesting  parameter, 

3(ln  C) _ 3(ln  C)  3 ( 1 n R) 

3(ln  L)  3(ln  R)  ' 3 ( 1 n L)  ’ 

can  indicate  the  sensitivity  of  cost  to  peak  demand.  The  term 
has  been  described  in  Table  5-1,  and  the  term  has  been  summa- 

rized in  Table  5-3.  Thus,  we  can  easily  derive  f ) by  the  multipli- 
cation  of  these  previous  coefficients.  Table  5-4  summarizes  these 
sensitivity  parameters. 
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Looking  at  k-mix  Type  4 and  k-mix  Type  1,  we  find  that  cost  sensi- 
tivity goes  up  5 times  as  the  capacity  of  generator  units  is  increased 
only  10  percent.  Thus,  it  appears  that  the  capacity  mix  is  very  cost 
sensitive  to  peak  demand. 

5.2.3  Interpretation  of  Cost  Equations:  Comparisons  between 

Models  with  Reliability  Variables  and  without 
Reliability  Variables 

2 

By  examining  Table  5-5,  the  logarithmic  model,  we  found  that  the  R 

2 

of  models  with  reliability  variables  are  slightly  higher  than  the  R of 

the  model  without  reliability  variables.  Thus,  the  explanatory  power  of 

the  model  is  improved  by  including  reliability  variables  into  the  cost 

model.  Note  that  the  most  meaningful  interpretation  of  the  model  with 

reliability  variables  is  to  assess  the  impact  of  reliability  on  cost, 

2 

even  though  the  improvement  of  R is  small.  The  statistical  hypothesis- 
testing showed  that  all  those  estimated  coefficients  are  statistically 
significant.  Looking  at  the  coefficients  of  total  electricity  consump- 
tion in  Table  5-5,  (i.e.,  , we  found  those  models  with  either 

LOLP  or  e(DNS)  have  higher  coefficients.  However,  for  K-mix  Type  4, 
such  coefficients  are  only  slightly  higher. 

2 

Looking  at  Table  5-6,  the  linear  model,  all  R are  above  99  per- 
cent. The  marginal  cost  of  electricity  production  for  the  model  without 
the  reliability  variable  also  moves  higher  as  each  generator  unit  is 

increased  by  10  percent,  that  is,  K-mix  Type  4.  Similar  to  the  log 

2 

model,  the  inclusion  of  reliability  in  the  model  does  not  improve  R 
substantially.  However,  it  does  enable  us  to  assess  the  impact  of 
reliability  on  cost.  By  pooling  the  data  from  all  four  mixes  (reported 


Table  5-5.  Comparisons  of  Log  Models  with 
and  without  the  Reliability  Variable 


Model  without 

Model  with 

Model  with 

Rel  iabil  it.y  Variable 

L0LP  Variable 

e(DNS)  Variab 

K- 

■mix  Type  1 

1 . 

Total  electricity 

consumption 

0.479 

0.5500 

0.6060 

2. 

Reliability  index 

- 

-0.0106 

-0.0150 

3. 

Constant  term 

•'  7.680 

6.8800 

6.1900 

4. 

R2 

99.67% 

99.89% 

99.82% 

K- 

■mix  Type  2 

1. 

Total  electricity 

consumpti on 

0.482 

0.5600 

0.5600 

2. 

Reliability  index 

- 

-0.0077 

-0.0086 

3. 

Constant  term 

7.640 

6.7500 

6.4600 

4. 

R2 

99.69% 

99.85% 

99.73% 

K- 

•mix  Type  3 

1. 

Total  electricity 

consumption 

0.535 

0.6200 

0.6800 

2. 

Reliability  index 

- 

-0.0140 

-0.0170 

3. 

Constant  term 

6.940 

5.9000 

5.1900 

4. 

R2 

99.49% 

99.82% 

99.65% 

K- 

•mix  Type  4 

1 . 

Total  electricity 

consumption 

0.450 

0.4800 

0.4600 

2. 

Reliability  index 

- 

-0.0340 

-0.0570 

3. 

Constant  term 

8.070 

7.7100 

8.7800 

4. 

R2 

99.55% 

99.97% 

99.97% 
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Table  5-6.  Comparisons  of  Linear  Models  with 
and  without  the  Reliability  Variable 


Model  without 

Model  with 

Model  with 

Reliability  Variable 

LOLP  Variable 

e(DNS)  Variab 

K- 

mix  Type  1 

1. 

Total  electricity 

consumption 

1 .66 

1.74 

1 .47 

2. 

Reliability  index 

- 

-715.60 

173.20 

3. 

Constant  term 

420,934 

404,700 

460,187 

4. 

P2 

99.89% 

99.91% 

99.97% 

K- 

mix  Type  2 

I. 

Total  electricity 

consumption 

1 .67 

1.56 

1 .53 

2. 

Rel iabi 1 ity  index 

- 

788.48 

171 .70 

3. 

Constant  term 

418,539 

440,745 

448,914 

4. 

R2 

99.90% 

99.93% 

99.99% 

K- 

mix  Type  3 

1 . 

Total  electricity 

consumption 

1.79 

1 .87 

1.45 

2. 

Reliability  index 

- 

-580.70 

293.86 

3. 

Constant  term 

360,662 

346,151 

428,984 

4. 

R2 

99.77% 

99.78% 

99 . 99% 

K- 

mix  Type  4 

1. 

Total  electricity 

consumption 

1.82 

1 .83 

1 .84 

2. 

Reliability  index 

- 

-184.70 

-21 .40 

3. 

Constant  term 

504,612 

500,796 

500,674 

4. 

R2 

99.99% 

99 . 99% 

99 . 99% 
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in  the  next  section),  we  find  that  the  reliability  variable  signifi- 
cantly improves  the  power  of  the  model. 

5.3  Comparison  of  Small  Variation  k-Mix  and  Large  Variation  k-Mix 

The  examination  of  specific  capital  mixes  involved  substantial 
multicollinearity.  Here,  we  analyze  the  results  when  the  data  are 
pooled.  In  addition,  in  order  to  explore  the  significance  of  k-mix 
variation  on  the  estimation  of  the  cost  equation,  we  define  another  set 
of  four  capacity  mixes  as  a large  variation  k-mix  group.  This  new  group 
contrasts  with  the  small,  variation  k-mix  group  described  in  Figure  4-3, 
whose  model  tests  were  reported  previously.  The  four  mixes  include: 

1.  Reference  mix  (Total  capacity  = 70  MW) 

Unit  B1  = 22  MW,  Unit  B2  = 30  MW,  Unit  A = 18  MW 

2.  Low  capacity  mix  (Total  capacity  = 60  MW) 

Unit  B:  = 21  MW,  Unit  B2  = 25  MW,  Unit  A = 14  MW 

3.  Lowest  capacity  mix  (Total  capacity  = 50  MW) 

Unit  Bx  = 20  MW,  Unit  B2  = 20  MW,  Unit  A = 10  MW 

4.  High  capacity  mix  (Total  capacity  = 80  MW) 

Unit  Bj  = 25  MW,  Unit  B2  = 35  MW,  Unit  A = 20  MW 

Given  the  problems  with  interpreting  results  for  specific  capital 

mixes,  we  proceed  directly  to  the  capital  mix  groupings.  The  results  of 

model  test  are  reported  in  Table  5-7.  First,  consider  models  with  and 

without  the  reliability  variable.  The  addition  of  L0LP  or  e(DNS)  en- 

? 

hance  the  explanatory  power.  The  reported  R are  substantially  in- 
creased by  including  reliability  in  the  model.  This  result  suggests 
that  parameter  cost  equations  reported  in  the  empirical  literature  could 
involve  biased  estimates. 

Comparing  the  small  variation  group  with  the  large  variation  group, 
we  find  the  estimated  intercepts  in  the  small  variation  group  are 


Table  5-7.  Comparisons  of  Parameter  Estimates  with 
and  without  Reliability  Variables 
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negative  when  reliability  is  included  as  an  explanatory  variable.  These 
negative  intercepts  are  unreasonable  from  the  standpoint  of  economic 
theory.  Thus,  it  may  be  inappropriate  to  pool  those  k-mixes  with  small 
variations  into  a single  group.  We  need  to  individually  estimate  each 
k-mix,  as  was  reported  in  the  previous  sections  of  this  chapter.  As  for 
the  large  variation  group,  the  intercept  coefficients  are  statistically 
significant  and  of  the  correct  sign. 

Turning  to  the  log  model  versus  the  linear  model  comparison,  we 
2 

note  that  all  R of  the  log  models  are  substantially  higher  than  those 
2 

R of  the  linear  models.  Thus,  the  log  model  has  greater  explanatory 
power,  perhaps  because,  as  Brown  and  Johnson  (1969)  note,  uncertainty  in 
demand  can  change  a constant  cost  technology  into  one  with  scale 
economies. 

Table  5-8  presents  the  reliability  elasticities  of  cost  estimated 
from  the  log  model.  The  three  estimated  coefficients  do  have  the  cor- 
rect signs.  For  the  large  variation  in  capacity  mix,  the  R2  of  the 
model  with  e(DNS)  as  the  index  is  highest.  As  to  the  reliability  elas- 
ticity of  cost,  cost  is  only  slightly  more  sensitive  to  the  e(DNS) 
reliability  index  than  to  LOLP.  Planners  need  to  integrate  into  their 
analysis  the  benefits  associated  with  small  reliability  improvements. 
The  results  reported  here  suggest  that  LOLP  and  e(DNS)  are  highly  cor- 
related for  the  particular  cost  and  demand  structures  used  here. 

5.4  Price  Setting  with  Reliability  Constraint 

In  this  concluding  section,  we  will  present  an  interesting  result 
of  this  simulation  test.  Suppose  a power  company  is  considering  whether 
to  adopt  time-of-day  (TOD)  pricing.  What  is  the  best  price  structure 
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for  the  community  subject  to  reaching  a given  level  of  reliability 
constraint?  Recall  that  the  capital  mix  is  assumed  to  be  fixed,  so  we 
are  considering  short-run  prices  (revenue  requirements  are  met  through 
lump  sum  fees).  Assume  the  electric  utility  sets  the  LOLP  reliability 
constraint  goal  to  be  less  than  2.5  days  per  year.  By  reviewing  Table 
5-9,  we  can  easily  locate  the  best  TOD  price  as  8.4<t  for  the  peak  and 
4.4<t  for  the  off-peak.  This  is  simply  because  the  benefit  is  maximized 
subject  to  the  reliability  constraint  at  this  particular  price  combina- 
tion. Of  course,  the  setting  of  optimal  reliability  constraint  involves 
the  estimation  of  cost  and  benefit  of  reliability.  The  estimation  of 
benefit  of  reliability  is  beyond  the  scope  of  this  study;  however,  this 
is  an  important  area  for  future  research. 
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CHAPTER  6 
CONCLUSIONS 


6.1  Summary 

This  study  has  arisen  out  of  a concern  on  the  part  of  the  public 
and  of  regulators  that  utilities  may  be  incurring  excessive  costs  for 
reliability.  On  the  other  hand,  these  groups  have  also  expressed  in- 
terest in  the  fact  that  currently,  utilities  are  cutting  deeply  into 
their  construction  programs  as  they  seek  to  curb  expenditures.  If 
present  trends  continue,  the  reliability  of  electric  services  could  be 
severely  jeopardized  in  the  early  1990s.  The  key  contribution  of  eco- 
nomic analysis  to  the  debate  is  that  the  cost  of  reliability  is  signifi- 
cant, requiring  decision  makers  to  choose  prices  and  capacity  so  as  to 
balance  additional  benefits  against  additional  costs  of  reliability. 

The  estimation  of  reliability  benefits  (avoided  outage  costs)  and 
supply  costs  is  necessary  for  a full  application  of  the  reliability 
optimization  model.  Instead  of  trying  to  capture  both  benefit  and  cost 
sides  of  alternative  reliability  levels,  this  study  sharpens  our  under- 
standing of  the  cost  of  reliability  improves  and  our  knowledge  about  the 
implications  of  alternative  reliability  indices  for  decision-making. 

This  study  differs  from  previous  economic  studies  in  three  res- 
pects. The  first  is  that  the  reliability  of  power  supply  becomes  an 
essential  variable  in  the  analysis  of  electricity  costs.  Two  relia- 
bility indices  widely  used  by  engineers,  LOLP  and  e(DNS),  are  analyzed 
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and  constructed.  The  second  contribution  of  this  study  is  that  demand 
and  supply  are  both  treated  as  stochastic  variables.  Most  previous 
studies  have  considered  "demand"  as  the  only  stochastic  element,  which 
limits  the  applicability  of  their  analytical  results.  The  third  aspect 
of  the  study  warranting  some  attention  is  its  use  of  pseudo  data.  The 
numbers  allow  us  to  explore  how  econometric  estimation  techniques  (and 
model  specification)  affect  estimates  of  the  cost  of  reliability. 

In  Chapter  3 we  provided  definitions  and  analyzed  implications  of 
alternative  reliability  indices.  Then  the  cost  characteristics  of  the 
utility  industry  were  reviewed  to  facilitate  the  development  of  the 
model.  Each  of  the  indices  was  shown  to  have  its  own  strengths  and 
weaknesses;  none  alone  provide  a complete  characterization  of  outage 
patterns  nor  of  the  value  of  lost  consumption.  The  estimated  indices 
are  only  significant  in  relative  rather  than  absolute  terms.  Since  the 
indices  are  generally  neither  comparable  nor  consistent,  they  yield 
different  rankings  of  systems  by  reliability. 

Due  to  lack  of  sufficient  and  reliable  data,  pseudo  data  based  on 
underlying  engineering  and  economic  considerations  were  used  in  this 
study.  In  Chapter  4,  we  presented  a simulation  model  to  generate  pseudo 
data.  We  then  used  alternative  regression  models  to  estimate  the  cost 
of  reliability  in  Chapter  5. 

The  model  we  proposed  is  a set  of  two  equations  for  reliability 
cost.  The  explanatory  variables  in  the  cost  equations  are  total  elec- 
tricity consumption  and  a reliability  index.  Peak  load  is  the  only 
explanatory  variable  in  the  reliability  equation.  We  separately  esti- 
mated both  linear  and  logarithmic  models  for  each  of  four  capacity 


106 


mixes.  The  two-stage  least-square  method  was  used  to  estimate  these 
equations. 

A new  concept,  reliability  elasticity  of  cost,  was  introduced  to 
measure  the  sensitivity  of  cost  to  reliability.  Those  estimated  elas- 
ticities shown  in  Table  5-1  are  all  much  less  than  unity.  Also,  analy- 
sis  of  standard  errors  of  estimated  coefficients  indicates  that  the 
e(DNS)  index  model  gives  the  best  estimate.  This  result  also  holds  for 
the  marginal  cost  of  electricity  production. 

6.2  Research  Recommendations 

As  discussed  previously,  the  estimation  of  reliability  benefits 
(avoided  outage  cost)  and  supply  cost  is  the  prerequisite  to  applying 
the  reliability  optimization  model.  However,  we  only  address  the  cost 
side  of  reliability  in  this  study.  As  to  the  benefit  side,  very  little 
hard  evidence  exists  in  economic  literature.  Thus,  the  problem  of 
reliability  benefit  is  a key  topic  for  future  research. 

In  addition,  load  management  techniques  can  be  used  to  assure  that 
peak  period  demand  is  within  the  capacity  limit.  Two  promising  innova- 
tive load  management  practices  are  peak  load  pricing  and  interruptible 
service  contracts.  The  former  has  been  throughly  examined  in  the  eco- 
nomic literature.  The  results  presented  here  show  that  higher  peak 
prices  increase  observed  reliability,  and  can  reduce  capacity  require- 
ments. The  latter  is  relatively  unexplored.  Thus,  it,  too,  is  worthy 
of  further  research.  The  essence  of  interruptible  service  pricing  is 
that  the  utility  has  the  option  of  curtailing  delivery  of  electricity  in 
a predetermined  order  when  available  supply  falls  short  of  demand.  The 
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customers  are  charged  different  prices  depending  on  their  position  in 
this  priority  ordering.  Also,  by  virtue  of  their  position,  they  con- 
front different  reliabilities  of  service,  and  the  influence  of  these 
reliabilities  on  demand  can  be  taken  into  account  by  system  planners. 

Finally,  one  worthwhile  point  for  future  research  is  that  public 
utility  commissions,  environmentalists,  and  various  citizen  groups  are 
asking  electric  utilities  to  build  smaller  generating  units  than  those 
built  in  the  past.  They  suggest  that  a higher  reliability  of  a 
portfolio"  of  smaller  units  would  result  in  a reduction  in  generating 
capacity  requirements,  which  could  result  in  investment  savings  that 
offset  the  loss  of  economies  of  scale  from  having  large  units.  However, 
others  argue  that  the  premise  of  economies  of  scale  is  still  alive, 
implying  that  large  units  are  more  economical  than  smaller  ones.  This 
controversial  issue  also  warrants  further  investigation.  The  results 
presented  here  are  offered  as  first  steps  in  more  appl ications-oriented 
analyses  of  utility  pricing  and  investment  policies. 
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